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1^0  ABSTRACT  (Continue  on  reverse  side  If  necessary  and  identity  by  block  number) 

\he  electrical  response  of  a copper-constantan  junction  to  shock  compres- 
sion has  been  studied  over  a pressure  range  from  145  kbars  to  360  kbars.  Four 
possible  sources  of  anomalous  response  were  found;  electrical  noise  due  to 
circuit  closure  at  impact,  shock  demagnetization  of  a ferromagnetic  material, 
high  local  temperature  at  the  junction  interface  due  to  shock  compression  of 
a surface  damage  layer,  and  two-dimensional  flow  in  the  pressurized  region  due 
to  pressure  relief  from  the  edges.  Using  a diffusion  welded  junction  the  emf  ^ 
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/^measured  is  within  20%  of  the  predicted  emf  based  on  normal  thermocouple 
/ response  to  shock  compression  temperatures.  In  a geometry  in  which  a guard 
ring  was  used  to  maintain  the  pressure  behind  the  shock  front,  the  observed 
voltage-time  profile  was  a step  with  constant  plateau.  In  a similar  configu- 
ration where  radial  pressure  relief  was  allowed  behind  the  shock  front,  the 
initial  step,  which  compares  directly  in  amplitude  to  the  above  results,  is 
followed  by  a marked  positive  ramping.  This  ramping  was  directly  correlated 
with  the  two-dimensional  flow  due  to  radial  pressure  relief.  The  results 
indicate  thay  any  fundamental  anomaly  which  exists  for  one-dimensional  strain 
is  less  than  10-20%  of  the  observed  signals  in  the  pressure  range  studied. 
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I INTRODUCTION 

Numerous  investigations  undertaken  since  1959  have  shown  that  when  a 

- shock  wave  passes  through  a junction  of  dissimilar  metals,  an  emf  is  observed 

I 1 

I j between  the  unshocked  extremities  of  the  metals.  Since  the  temperature  is 

I ^ 

f ; elevated  in  the  shocked  region,  an  emf  is  expected  if  the  metal  pair  behaves 

I as  a normal  thermocouple.  In  every  case  the  observed  emf  was  of  the  appro- 

\ ■ 

priate  sign  but  greater  in  magnitude  than  would  be  expected  from  normal 
thermoelectric  response  to  shock  compression  temperatures.  In  most  studies 
the  emf  was  observed  to  be  between  3 and  6 times  too  large.  A great  deal  of 
confusion  exists  as  to  the  source  of  the  abnormally  large  response  and  the 
cause  of  the  general  nonreproducibility  among  experiments. 

I Initially  such  studies  were  motivated  by  the  need  for  a transducer 

f 

I which  would  continuously  indicate  shock  strength  as  a function  of  time. 

I Attempts  were  made  to  calibrate  the  voltage  output  with  pressure  or  temper- 

! ature  but  nonreproducibility  among  investigators  has  made  this  impossible. 

I In  some  investigations  primary  emphasis  was  placed  on  determining  the  cause 

; of  the  abnormally  large  emf.  Some  progress  has  been  made  but  no  general  reso- 

i. 

i lution  of  the  problem  now  exists. 

The  objective  of  this  study  was  to  bring  some  order  into  the  study  of 
this  effect  by  carefully  measuring  the  response  of  a bimetallic  junction  to  a 
step  pressure  pulse  using  a light  gas  gun  and  modern  high  resolution  recording 
equipment.  The  experimental  program  developed  in  two  phases,  the  first  of  which 
was  to  develop  an  experimental  configuration  in  which  conditions  such  as  junc- 
tion quality  and  shock  pressure  could  be  completely  determined.  It  was  also 


L ^ 
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necessary  to  show  that  no  parasitic  signals  were  inadvertently  being  measured 
due  to  sources  othar  than  the  junction  response.  The  objective  of  the  second 
phase  of  this  program  was  to  use  the  experimental  configuration  so  developed 
to  study  the  junction  response  over  a range  in  pressures  from  145  kbars  to 
360  kbars,  and  to  compare  it  to  anticipated  response  based  on  conventional 
static  calibrations. 

Throughout  the  course  of  this  investigation,  the  primary  question  to 
be  answered  in  light  of  the  previous  state  of  confusion  was  whether  the 
anomalously  high  emf  and  nonreproducibility  previously  reported  were  due  to 
fundamental  principles  or  experimental  artifacts.  The  results  reported  here 
indicate  not  only  that  anomalously  large  emfs  were  indeed  due  to  experimental 
artifacts,  but  also  that  the  simple  measurement  of  shock  compression  temper- 
ature with  a thermoelectric  circuit  may  be  possible. 
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CHAPTER  1 

THE  SEEBECK  EFFECT  AND  APPLICATION  TO  THE  SHOCK 
COMPRESSED  BIMETALLIC  JUNCTION 

1.1  Thermoelectric  Power  of  Metals. 

The  present  theory  of  thermoelectricity  is  neither  exact  nor  complete. 
The  basic  relationship  between  causes  and  effects  involved  are  embodied  in 
the  simple  linear  Onsager  relations,  but  an  exact  theoretical  understanding 
of  the  coefficients  is  precluded  by  complexities  of  real  materials  on  the 
atomic  scale.  The  thermoelectric  phenomena,  which  include  the  Seebeck, 

Peltier  and  Thomson  effects,  are  combinations  of  more  than  one  process  invol- 
ving electron  diffusion  as  well  as  electron-electron,  electron-phonon,  and 
electron-impurity  interactions.  Although  each  of  these  processes  can  be 
treated  in  a general  way  with  the  use  of  modern  quantum  mechanics  and  statis- 
tical mechanics,  it  is  not  possible,  except  in  rare  cases,  to  predict  quanti- 
tative results  as  accurately  as  they  may  be  measured.  At  this  point  in  the 
development  of  the  general  theory,  experimental  results  of  thermoelectric 
phenomena  are  being  interpreted  in  terms  of  the  general  formalism  in  order  to 
gain  a fuller  understanding  of  both  the  theory  and  the  specific  properties  of 
the  materials.  The  objective  of  the  following  treatment  is  to  give  a working 
understanding  of  mechanisms  known  to  be  active  in  the  materials  of  interest 
in  order  to  evaluate  the  experimental  results  of  the  current  investigation  in 
terms  of  results  of  previous  work  under  normal  laboratory  conditions.  This 


discussion  will  be  limited  to  the  Seebeck  effect  in  metals  with  particular 
attention  given  to  constantan,  an  alloy  of  55%  copper  and  45%  nickel.* 

The  Seebeck  effect  or  Seebeck  potential  is  the  electric  potential  dif- 
ference, A<1),  which  exists  between  the  ends  of  a material  in  which  a temper- 
ature gradient  is  maintained.  The  common  measurement  of  this  effect  is  made 
using  a circuit  as  in  Fig.  1.1  in  which  two  materials  with  different  proper- 
ties are  used  and  no  current  is  allowed  to  flow.  The  thermoelectric  power  of 
the  couple  is  defined  by 

S.  R = lim  (A({»/AT).  (1.1) 

It  is  not  possible  under  normal  conditions  to  measure  the  potential  difference 
due  to  a temperature  gradient  in  a single  material  because  the  measurement 
itself  involves  other  materials.  It  is  possible,  however,  to  define  an 
absolute  thermopower  by  the  relation 


E = (1.2a) 

where  E is  the  electric  field  in  the  material  and  vT  is  the  temperature 
gradient.  In  general  ^ is  a second  rank  tensor.  In  any  further  treatment 
here,  except  where  noted,  it  will  be  assumed  that  the  materials  of  interest 
are  isotropic,  and  S can  be  written  as  a scalar.  Since 


Eq.  (1.2a)  becomes 


-V4) 


d(f)  = -S  dT 


where  (}>  is  the  electric  potential. 


(1.2b) 


*The  name  constantan  is  commonly  used  for  a class  of  copper-nickel 
alloys  ranging  in  composition  from  60%-40%  to  50%-50%,  and  sometimes  inclu- 
ding 1%  manganese. 


Fig.  1 .1 .--Thermoelectric  Circuit  through  Temperatures 
T and  T + AT. 
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The  absolute  thermopower  and  associated  Seebeck  effect  represent  only 
one  part  of  a more  complex  interdependence  among  temperature  gradients,  elec- 
tric potential  gradients,  heat  flux,  and  electric  current.  The  formal 
relationships  among  these  quantities  are  known  as  the  Onsager  relations  and 
will  be  treated  here  within  the  framework  of  irreversible  thermodynamics  (see 
Callen^),  In  the  theory  of  irreversible  thermodynamics  it  is  assumed  that  the 
intensive  parameters  may  be  defined  as  they  are  in  equilibrium  thermodynamics 
as  the  appropriate  derivatives  of  the  extensive  parameters,  and  that  there 
exist  generalized  forces  called  affinities  and  conjugate  fluxes. 

In  a system  in  which  no  volume  changes  are  allowed,  one  can  write  a 
general  expression  for  the  differential  of  the  local  entropy  density  as 

dC  = I du  - ^ dn  - I (^)dn^  (1.3) 

where  C is  the  local  entropy  density,  u is  the  local  energy  density,  p is 
the  electrochemical  potential  per  electron,  and  n is  the  number  of  electrons 
per  unit  volume.  The  summation  term  is  included  to  take  into  account  any 
other  components  such  as  atomic  nuclei  which  make  up  the  solid. 

From  Eq.  (1.3)  one  can  define  an  entropy  flux  density  as 


- — — "l 

T '^u  ■ T '^n’ 


(1.4) 


where  and  are  the  energy  and  electron  flux  density,  respectively. 

The  other  components,  nj^,  are  neglected  in  the  further  development  as  they 
are  considered  inmobile  and  thus  do  not  contribute  fluxes  to  Eq.  (1.4).  In 
this  way  we  are  neglecting  effects  due  to  material  motion  in  this  treatment. 
The  rate  of  local  entropy  production  is  equal  to  the  amount  of  entropy  leaving 
the  region  plus  the  rate  of  local  entropy  increase.  Therefore, 
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1 
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k = ^ + v.j 
^ 3t 


(1.5) 


where  t is  the  rate  of  local  entropy  production,  and  d^/dt  is  the  rate  of 
local  entropy  increase.  The  appropriate  continuity  equations  for  energy  and 
number  of  electrons  are 


(1.6a) 


0 * . 


(1.6b) 


From  Eq.  (1.3)  the  first  term  of  Eq.  (1.5)  can  be  evaluated  as 

^ - i ly.  li  In 

3t  " T at  ■ T 3t  ' 


(1.7) 


By  taking  the  divergence  of  Eq.  (1.4)  and  using  (1.7),  Eq.  (1.5)  becomes 


+ (V  (T))*Jy  - (V  (T))-Jn 


(1.8) 


The  first  two  terms  are  zero  from  the  continuity  Eqs.  (1.6),  so  that 


C * (V  (|))*Ju  - (V  (f))-Jn  • 


(1.9) 


In  a continuous  system  the  affinities  of  irreversible  thermodynamics 
are  defined  as  gradients  of  the  entropy  representation  intensive  parameters. 
In  this  system  the  affinities,  or  generalized  forces,  are  V(l/T)  and 
V(y/T).  The  fluxes  are  and  -J^.  The  preceding  argument  is  given  to 
identify  the  form  of  the  fluxes  and  affinities  so  that  the  general  relation 
between  fluxes  and  affinities  may  be  introduced.  If  it  is  assumed  that  the 
fluxes  at  a given  instant  depend  only  on  the  affinities  at  that  instant,  the 
fluxes  can  be  expanded  in  a power  series  of  the  instantaneous  affinities  as 


’i  " I ^i/j  2!  ‘■ijk''j^k  + • • • 


(1.10) 


where  F.  represents  the  affinity  conjugate  to  the  flux  . Taking  only 

J J 


the  linear  term,  Eq.  (1.10)  becomes  for  the  present  case 


(1.11a) 

'’u  ° '■22''4'  • 

(1.11b) 

It  is  easier  to  work  with  the  flow  of  heat  than  with  the  flux  of  total 
internal  energy.  Following  the  same  procedure  used  to  define  the  fluxes  in 
Eq.  (1.4),  a heat  flux  can  be  defined  as 

Jq  = TJ^  , (1.12) 

and  by  Eq.  (1.4), 


The  rate  of  entropy  production  can  now  be  written  as 

i = V(|)-Jp  - (^)VyJ^  . (1.14) 

The  fluxes  and  affinities  are  now  reselected  and  Eq.  (1.11)  becomes 


■« 

M 

C 

1 

L^^(^)Vy  + 

(1.15a) 

II 

L2i(})Vy  + L22'^(t)  • 

(1.15b) 

The  Onsager  theorem,  where  no  externally  applied  magnetic  field  is 

present,  states  that  L12  = L2^.  For  a discussion  of  the  effects  of  an  exter- 

2 

nal  magnetic  field  on  this  relation  see  Ziman,  Chapter  XII. 

The  coefficients  L . . may  be  determined  in  terms  of  familiar  quanti- 

* J 

ties  such  as  heat  conductivity  by  comparing  special  cases  of  Eq.  (1.15)  with 


9 

established  physical  laws.  Using  the  Onsager  theorem,  one  can  then  discuss 

the  more  complex  combinations  of  fluxes  and  forces  in  terms  of  well  known 

phenomena.  As  an  example,  one  can  find  a relation  between  the  coefficients 

and  the  heat  conductivity  k by  solving  Eqs.  (1.15)  under  the  conditions 

that  the  electric  current  J is  zero  and 

n 

K = -Jp/VT  . (1.16) 

The  resulting  relation  is 

tc  = (L^^L22  - • (1-17) 

In  this  way,  the  three  independent  coefficients  can  be  determined  in  terms  of 
the  electric  conductivity  o,  the  heat  conductivity  k,  and  the  absolute 
thermopower  S.  Equation  (1.15)  can  now  be  written 

e 

Jq  = -(^)(})Vu  + (tW  + T2K)V(i)  , (1.18b) 

where  e (<0)  is  the  charge  of  an  electron. 

The  thermopower  expression,  Eq.  (1.2b),  can  be  regained  from  Eq.  (1.18) 
by  using  the  first  of  the  two  expressions  and  letting  the  current  be 
zero,  which  yields 

0 * ^|)dy  - (I^)d(|)  (1.19) 

or 


-S  dT  . 


(1.20) 
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The  electrochemical  potential,  y,  can  be  written  as  a sum  of  two  components, 
the  electric  part  and  the  chemical  part,  as 

y = ^c  ’ (1-21) 

where  y^  = e4),  and  y^  depends  on  temperature  and  electron  concentration. 
Equation  (1.20)  becomes 

(i)  dy^  + del)  = -S  dT  (1.22a) 


or 


d(l)  = -S  dT  - (1)  dy^  . (1.22b) 

The  extra  term  in  Eq.  (1.22)  as  compared  with  Eq.  (1.2b)  gives  rise  to  the 
contact  potential  or  the  difference  in  electric  potential  across  a material 
boundary.  This  can  be  seen  by  integrating  Eq.  (1.22b)  along  a small  path 
across  the  boundary  between  two  metals  I and  II  so  that 

♦ll  -♦!  = -“J)  • ''-23) 


Since  temperature  is  continuous,  the  other  term  vanishes  in  the  limit  as  the 
integration  path  length  goes  to  zero.  This  contact  potential  is  not  measured 
explicitly  in  any  thermoelectric  circuit  as  can  be  seen  in  the  following 
example.  Consider  the  circuit  of  Fig.  1.1.  The  materials  A and  B have  abso- 
lute thermopowers  and  Sg,  respectively.  The  difference  in  the  electro- 
chemical potential  between  point  a and  b is  given  by  the  integration  of 


Eq.  (1.20). 


a b 
y - y 


= e 


(-S^)dT'  + e 


T+AT  'a 

(-Sg)dT'  + e 


(-S^)dT* 


(1.24) 


T+AT 


Since  the  temperature  and  the  material  are  the  same  at  a and  b,  the  chemical 
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potential  u is  the  same  at  a and  b,  so  that  from  Eq.  (1.21) 

a b a b 

y - y * Pg  - Pg  = e(<))  - (J>  ) 


(1.25) 


and  therefore 


T+AT 


(-SB)dT'  + 


(-S;^)<1T' 


(1.26) 


T+AT 


This  can  be  compared  with  the  defining  equation  for  the  thermoelectric  power 
of  such  a couple,  (Eq.  1.1).  By  combining  the  integrals,  Eq.  (1.26)  becomes 


T+AT 


A(|)  ■ (j)®  - ({i^ 


(S,  - Sj,)dT 


(1.27) 


and  by  comparison  with  Eq.  (1.1) 

Vb  “ " ^B  • 

It  is  important  to  note  here  that  the  current  which  flows  in  a thermo- 
electric circuit  is  in  the  opposite  direction  from  that  of  the  net  external 
electric  field.  This  has  led  to  a great  deal  of  confusion  as  to  the  sign  of 
the  thermoelectric  voltage.  As  is  the  case  in  any  emf  source  such  as  a vol- 
taic cell  or  thermocouple,  there  is  an  opposing  force  within  the  source  which 
drives  the  current  in  the  opposite  direction  from  that  of  the  external  elec- 
tric field.  In  a thermocouple  the  confusion  seems  to  arise  because,  unlike 
other  emf  sources,  it  is  not  obvious  that  the  source  is  anything  but  a resis- 
tive element  of  a circuit.  In  a circuit  such  as  in  Fig.  1.1  the  entire 
circuit,  apart  from  the  measuring  apparatus,  can  be  considered  an  emf  source 
with  an  associated  internal  resistance  much  like  a battery. 


A 


The  two  basic  mechanisms  which  give  rise  to  the  absolute  thermopower 
in  metals,  and  are  represented  by  the  Onsager  relations,  are  thermal  diffusion 
of  electrons  and  what  is  called  phonon  drag.  The  thermopower  is  usually 
written  as  a sum  of  two  terms 


L 


s = S . + 

d g 


(1.29) 


where  and  refer  to  the  diffusion  and  phonon  drag  thermopowers, 
respectively.  The  diffusion  process  accompanies  heat  transport  by  the  mobile 
electrons,  and  the  phonon  drag  term  is  due  to  an  interaction  of  the  electrons 
with  heat  transport  through  the  crystalline  lattice. 

3 

The  phonon  drag  thermopower  is  negligible  above  room  temperature. 

The  diffusion  part  of  the  thermopower  may  be  discussed  in  terms  of  Boltzman 


transport  theory.  For  a more  complete  discussion  of  this  subject  see  Ziman 
4 


or  Mott  and  Jones.  The  significant  result  is  that  the  parameters  o,  k,  and 
S can  be  determined  in  terms  of  the  Fermi  energy,  the  Boltzman  constant, 
k,  and  the  temperature,  T.  It  can  be  shown  that 


a = e K_ 


K = (1/T)(K2  - K^K'^) 


(1.30a) 

(1.30b) 

(1.30c) 


where 


((e-c)"  2i|i  + ^ kV  4-  [(£-0"  2^  * 


. .} 


3e 


(1.31) 


Here  a(e)  is  to  mean  the  electrical  conductivity  which  one  would  calculate 
if  the  Fermi  energy  of  the  metal  were  e.  The  Mott  expression  for  the  thermo- 
power is  obtained  by  substituting  K.j  and  into  Eq.  (1.30c)  and 


r 


retaining  only  the  terms  to  first  order  in  kT/;.  The  result  is 


.2.2 


. _ ir^k^T  ,81na(e)x 

~ 3e  ' ■ 


3e  'e=^ 


(1.32) 


It  has  been  shown  that  for  a free-electron  metal  the  logarithmic  derivative 
in  Eq.  (1.32)  is  just  1/^  for  T«9p  and  3/c  for  T > 6^  where  9p  is 
the  Debye  temperature.  Equation  (1.32)  becomes 


ir^k^T 


3es 


T«9r 


(1.33) 


■n\h 


ec 


T > 9. 
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(1.34) 


This  expression  gives  a thermopower  of  about  1 Mvolt/°K  at  room  temperature 
for  reasonable  values  of  the  Fermi  energy.  The  linear  dependence  of  S on 
T has  been  shown  for  high  temperatures  (above  9^)  in  most  metals.  Although 
more  accurate  expressions  for  the  thermopower  can  be  derived  by  considering 
the  electron-lattice  interaction  more  carefully,  one  cannot  at  this  time  pre- 
dict the  thermopower  quantitatively  to  the  accuracy  of  experimental  results 
above  room  temperature. 


1.2  Shock  Compression  in  a Thermoelectric  Circuit. 

Consider  the  experimental  configuration  illustrated  in  Fig., 1.2.  A 
shock  wave  is  propagating  away  from  the  junction  in  both  copper  and  con- 
stantan.  It  will  be  assumed  that  no  pressure  relief  takes  place  from  the 
edges  behind  the  shock  although  no  means  of  support  is  shown.  The  junction 
will  be  considered  as  a perfect  plane  dividing  the  two  materials.  The  fol- 
lowing discussion  will  be  based  on  this  rather  unrealistic  configuration  and 
a comparison  will  be  made  later  with  the  experiments  actually  performed. 


I 
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Junction 


Copper  J,  Constonton 


A B 


Fig.  1 .2. --Idealized  Configuration  for 

Bimetallic  Junction  Experiments. 

The  shock  waves  are  produced  by 
plane  impact  of  the  two  components. 


i 


The  measured  quantities  in  this  idealized  experiment  are  the  emf 
between  A and  B (Fig.  1.2),  and  the  velocity  of  the  impactor.  The  velocity 
of  the  impactor  is  used  with  the  Rankine-Hugoniot  jump  conditions,  and  known 
empirical  relations  between  shock  velocity  and  particle  velocity  for  each 
material,  to  obtain  the  pressure  behind  the  shock  front.  Pressures  and 
temperatures  obtained  for  various  impactor  velocities  will  be  discussed  in 
the  following  section.  At  this  point  it  will  be  assumed  that  the  temperature 
is  elevated  in  the  metals  due  to  shock  compression  so  that  an  emf  of  thermo- 
electric origin  is  expected. 

The  predicted  thermoelectric  emf  can  be  obtained  from  an  integration 
of  Eq.  (1.2b),  where  S = S(T,P)  is  a function  of  temperature  and  pressure. 
Figure  1.3  illustrates  the  temperature  profile  for  this  idealized  experiment. 
Since  S(T,P)  will  be  different  for  the  two  metals,  the  integration  from  A 
to  B of  Eq.  (1.2b)  can  be  written  as 


‘*’B-A  - “^B  ■ '*’A 


-jWT'.'’) 

^0  . , 
region  I 


dT'  - 


Sc„(T'.P)  dT' 
* region  II 


Sc,(T-,P)  dT- 


(1.35) 


region  III 


where  T^  is  the  temperature  of  the  materials  before  the  experiments,  and 
Tj  is  the  temperature  at  the  junction  measured  in  degrees  Kelvin.  The  sub- 
scripts Cu  and  Co  refer  to  copper  and  constantan,  respectively,  and 
is  defined  as  the  potential  difference  (Jig  - (J)^.  If  S(T,P)  were  known  for 
each  region  along  with  the  functional  relationship  between  pressure  and 


Temp 
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Fig.  1 .3.--Tenperature  Profile  for  Idealized  Experiment. 
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temperature,  the  integration  could  be  performed  and  the  results  compared 
directly  with  experimental  results.  However,  empirical  relations  which  are 
available  for  thermopower  as  a function  of  pressure  and  temperature  involve 
particular  thermoelectric  circuits,  so  that  S(T,P)  can  be  obtained  only  by 
careful  examination  of  these  circuits. 

To  understand  the  dependence  of  the  thermopower  on  pressure,  consider 
the  thermoelectric  circuit  of  Fig.  1.4  made  up  of  a wire  of  a single  material. 
Using  Eq.  (1.2b)  the  measured  emf  can  be  written  as 

<Pc_Q  = - S(T',P)  dT'  (1.36a) 


1 2 0 

- S(T',P=0)dT'  - |s(T’,P=P^)dT'  - S(T' ,P=0)dT' . (1 

T-  T,  T. 


Combining  terms  the  equation  becomes 


(1.37) 


(1.38) 


Fia.  1 .4. --Thermoelectric  Circuit  through  Pressurized  Region. 


"^C-D  ^ -nP-)*  (T2-T1 ) 


(1.39b) 


5 6 

From  experiments  done  by  Bridgman  up  to  12  kbars  and  Bundy  up  to  72  kbars, 
the  constant  n can  be  evaluated  for  both  copper  and  constantan  as 

Tico  “ 0.026  yv/kbar°K 
= -0.002  yv/kbar°K  . 

Equation  (1.35)  can  now  be  simplified  in  terms  of  this  pressure  dependence  as 


Since  in  the  shock  front  the  system  is  not  in  thermodynamic  equilibrium,  it 
is  not  clear  that  a temperature-pressure  dependence  is  meaningful.  However, 
as  an  approximation,  the  irreversible  thermodynamic  definition  of  temperature 
as  the  local  derivative  of  internal  energy  will  be  assumed,  and  a linear 
dependence  of  temperature  on  pressure  will  be  used  in  the  form 


T - T. 


Tp-'^o 


(1.41) 


Here  Tp  and  Pp  are  the  equilibrium  values  behind  the  shock.  Equation 
(1.40)  becomes 


^Cu^Cu 


(T'-TJ 


'Cu''o' 


'CuPcu'^^’ 


. ’^Co'^Co  (Tf.  -TJ  - J \o'’co‘^^' 


'Co  '0 


KJ 

- [S(.y(T*,0)  - S(.^(T',0)]dT' 


(1.42) 
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where  = Pj,^  = final  shock  pressure;  and  are  the  shock  temper- 

atures in  copper  and  constantan  respectively.  The  first  and  third  terms  are 
integrations  over  the  shock  front  where  the  above  temperature-pressure  depen- 
dence applies.  The  second  and  fourth  integrals  are  along  the  temperature 
gradient  due  to  the  difference  in  shock  compression  temperatures  of  the  two 
metals  where  the  pressure  is  constant.  This  expression  can  now  be  simplified 
as 


4* 


B-A 


■\u^Cu^^J  ■ 


Cu  0\  4.  „ n It  Co  "^Ox 
2 ■ 2 ^ ’^Co’^Co'^J  ■ 2 " 2 ^ 


(1.43) 


'E-F 


/ 


(1.44) 
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where  Scu_co(^'»0)  = thermopower  defined  for  a 

copper-constantan  thermocouple  at  zero  pressure  as  in  Eq.  (1.1).  The  cor- 
rection calculated  in  this  way  for  the  pressure  effect  corresponds  to  about 
an  8%  decrease  in  the  expected  emf  at  300  kbars.  It  remains  only  to  find  the 
temperatures  T^,^,  Tj.^,  and  Tj  as  functions  of  other  shock  parameters  and 
Scy-co(T)  in  terms  of  a known  empirical  relation  to  estimate  the  emf  observ- 
able in  shock  experiments.  The  temperatures  involved  will  be  treated  sepa- 
rately in  section  1.3  and  ^^.^  ^.^(T)  can  be  easily  obtained  by  the  following 
argument. 

Consider  the  thermoelectric  circuit  of  Fig.  1.5.  The  emf  at  the 
potentiometer  is  given  by 


from  Eq.  (1.1)  so  that 
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T 

^Cu-Co^^^  ^ ^ ^ 

^0 

It  is  assumed  here  that  is  a constant,  which  implies  that 
depends  on  the  reference  temperature.  The  materials  used  in  these  experiments 
were  tested  for  thermopower  in  this  type  of  circuit  with  = 23  ± 2°C  which 
was  the  same  as  for  the  shock  experiments. 

The  emf  was  measured  for  the  materials  used  here  by  a method 

similar  to  that  of  Nagy  and  Toth.^  An  empirical  relation  for  (|)(T)  was 
obtained  in  the  form 

‘*’E-f(‘^“'^o^  " Q(T-Tq)  + R(T-Tq)^  (1-46) 

where  Q and  R are  constants.  The  data  and  the  values  obtained  for  Q 
and  R are  given  in  Appendix  A.  The  thermoelectric  power  of  this  couple  can 
now  be  written  as 

Scu-Co<T>  ' • (1.47) 

In  summary,  the  calculated  emf  for  this  idealized  shock  experiment 
can  be  obtained  by  evaluating  Eq.  (1.43)  where  and  can  be 

determined  in  independent  experiments  and  T^,  Tj,  T^.^,  T^.^,  and 
are  known  or  can  be  estimated  from  measured  quantities.  A method  for  esti- 
mating Tq,  Tj,  and  Tj.^  is  given  in  the  next  section.  This  idealized 
experiment  will  be  compared  to  the  actual  experiments  when  they  are  described 
in  Chapter  4. 
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1.3  Temperature  in  Shocked  Copper-Constantan  Couple. 

In  this  section  an  estimate  will  be  made  of  the  temperatures  T^.^, 
Tcu»  and  Tj  indicated  in  Fig.  1.3,  based  on  a modified  Christian  and  Walsh 

O 

technique.  The  temperature  at  the  junction  will  be  treated  in  terms  of  an 
initial  step  temperature  profile  with  subsequent  heat  flow  to  produce  a fixed 
temperature  at  the  interface  which  depends  on  the  diffusivities  of  the 
materials. 

The  temperature  rise  due  to  shock  compression  of  a material  may  be 
estimated  by  inverting  and  integrating  the  expression 

d£  - (||),dT  (1.48) 

or 

dC  = dT  + ic^  Ydv  (1.49) 

where  C is  the  entropy  density  and  y is  the  Gruneisen  parameter. 

In  a shock  process  there  is  entropy  production  due  to  plastic  defor- 
mation as  well  as  hydrodynamic  compression,  so  that 

dC  = d?^  + dC2  (1.50) 

where 

d^i  = Y T (dv  - V ^)  (plastic  deformation^)  (1.51) 

and 

(^x”^x  ^ V - V dP 

^ 2 Cl  - (-  -d^)]dv  (hydrodynamic)  (1.52) 

"hel 

where  p^  is  the  component  of  pressure  in  the  shock  propagation  direction 
and  the  subscript  HEL  refers  to  the  Hugoniot  elastic  limit.  In  the  first 


expression  t is  the  maximum  resolved  shear  stress,  and  y is  the  shear 
modulus.  The  second  expression  is  the  direct  result  of  the  Rankine-Hugoniot 
relation 


24 


(E-E„)  . (P,  * P„)(v„  - v)/2 


0.53) 


and  the  general  thermodynamic  relation 


dC  = I dE  + Y dv 


(1.54) 


Using  the  expression  for  dE  along  with  Eq.  (1.49)  and  differentiating  with 
respect  to  v,  one  obtains 


dv  V ' 3 Cy  u dv^ 


* <'’x  - '’x 


*HEL 
P - P 


(-  Tlir)] 


(1.55) 


To  integrate  this  differential  equation  one  must  assume  a dependence  of  y 
and  Cy  on  volume.  For  this  estimation  it  will  be  assumed  that  y/v  and 
Cy  remain  constant.  The  constant  y/v  will  be  designated  as  G,  thus 


6 = y/v  = y^/v^  . 


(1.56) 


The  right  side  of  Eq.  (1.55)  will  be  designated  as  F(v)  so  that  the  equa- 
tion can  be  written  simply  as 


|i*GT  = F(v) 


(1.57) 


This  equation  has  the  solution 


T = t e-®’  I e^^'pCvOdv'  . 


(1.58) 


Integration  of  this  expression  was  done  numerically  by  the  application  of 
Simpson's  Rule.  The  pressure  as  a function  of  volume  along  the  Hugoniot  is 
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calculated  point  by  point  using  published  U^-Up  relations  and  the  Rankine- 
Hugoniot  jump  conditions.^®  (See  Appendix  B.) 

Work  hardening  during  the  shock  process  could  be  included  by  letting 
T be  a function  of  volume.  If  one  assumes  the  same  dependence  of  t as 
would  be  observed  for  quasi  static  strain  to  low  pressure,  t has  the  form 

T = Tq  + j Cln(v/VQ)]  (1.59) 

which  is  derived  from  the  single  yield  equation 

Y = Yq  + H(eg  - e^^).  (See  Appendix  B)  (1.60) 

Here  is  strain  in  uniaxial  stress,  Y is  the  yield  stress,  and  H is  a 
constant.  In  Appendix  B a comparison  is  made  of  temperature  estimates  with 
and  without  a work  hardening  correction.  The  correction  is  only  about  1%  at 
300  kbars.  The  values  of  H measured  here  are  15.5  and  3.74  for  copper  and 
constantan,  respectively.  The  Hugoniot  elastic  limit  for  copper  was  taken 
from  published  measurements  of  elastic  precursor  amplitudes.  In  constantan, 
the  Hugoniot  elastic  limit  was  estimated  from  the  quasistatic  yield  stress 
and  elastic  constants  obtained  from  ultrasonic  sound  speed  measurements.  For 
a linear  elastic  isotropic  solid 

(Ref.  9)  (1.61) 

Xhel 


where  Y is  again  the  yield  stress  in  uniaxial  stress  and  v is  Poisson's 
ratio.  Appendix  B includes  the  values  for  all  of  the  parameters  necessary 
in  these  calculations,  as  well  as  resultant  temperature  estimates  for  various 
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pressures.  The  temperature  estimates  for  shocked  copper  of  Rice,  et  a1 
are  also  listed  in  Appendix  B for  comparison. 

In  the  region  very  near  the  junction,  a step  temperature  profile 
exists  just  after  the  shock  transit  through  the  junction  as  illustrated  in 
Fig.  1.6.  The  temperature  at  the  interface  is  time  independent  and  given  by 
Carslaw  and  Jaeger^ ^ as 


T 


J 


^Cu  - ^Co 
1 + B 


(1.62) 


where  B is  given  by 


B 


/^J/2  ^ 


f^xl/2 

■^Co 


(1.63) 


In  this  expression  k and  D are  heat  conductivity  and  diffusivity,  respec- 
tively. The  approximation  holds  where  the  density  and  specific  heat  are 
nearly  equal,  which  is  true  for  copper  and  constantan,.  The  conductivities 
for  copper  and  constantan  at  100°C  are  3.88  watts/cm°C  and  0.212  watts/cm°C, 
respectively,  so  that  B,  which  will  be  assumed  independent  of  temperature  is 


B = 4.3 


(1.64) 


The  temperature  at  the  junction  is  now 


T 


J 


(1.65) 


The  temperatures  necessary  to  predict  the  thermoelectric  response  of 
a shocked  junction  of  copper  and  constantan  are  estimated  by  using  Eq.  (1.58) 
for  the  temperature  behind  the  shock,  Tj.^  and  T^.^,  and  by  using  Eq.  (1.65) 
for  the  junction  temperature  Tj. 


27 


r 


I 

I 

I 


Tcu  I 


I 

I 


x=o 


Fig.  1 .6.— Temperature  Profile  Near  Junction  Just  after  Shock 
Compression.  The  dotted  line  represents  the  tem- 
perature profile  due  to  heat  conduction.  The 
temperature  Tj  is  independent  of  time. 
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CHAPTER  2 

REVIEW  OF  PREVIOUS  WORK 

Experimental  results  from  previously  reported  studies  of  junction 
response  to  shock  compression  are  summarized  in  Table  2.1.  The  materials 
used  and  the  method  of  forming  the  junction  are  listed,  as  well  as  the  mea- 
sured emf  and  shock  pressure.  Entries  are  listed  in  chronological  order  by 
reporting  date.  In  some  cases  the  authors  have  compared  the  measured  emf  to 
a calculated  emf  based  on  a normal  thermoelectric  response  to  shock  loading. 
This  calculated  value  is  listed  for  comparison  as  it  appeared  in  the  original 
reports.  In  some  cases  where  the  thermopower  is  known  and  this  comparison 
was  not  made,  we  have  included  a value  in  parentheses  for  comparison.  Since 
objectives  and  experimental  procedures  vary  widely  among  investigations,  a 

few  comments  on  each  study  are  given  here. 

1 2 

Jacquesson  reported  voltage-time  profiles  in  experiments  designed 
to  measure  the  pressure  profile  in  copper  by  inserting  a thin  constantan 
ribbon  between  pieces  of  copper.  He  concluded  that  the  emf  was  very  large 
compared  with  normal  thermoelectric  response,  but  that  some  indication  of 
the  pressure-time  profile  could  be  obtained  qualitatively.  The  actual  emf 
trace,  however,  includes  large  positive  and  negative  excursions  and  in 

general  cannot  be  correlated  with  the  expected  pressure-time  profile. 

13 

Ilyukin  and  Kologrivov  conducted  a series  of  experiments  in  which 


a disk  of  nickel  was  soldered  to  the  back  of  a larger  disk  of  copper  with  no 
radial  support.  The  emf  values  observed  correspond  to  a temperature  of 


TABLE  2.1.— Results  of  Previous  Investigations 
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1000°C  for  a set  of  experiments  in  which  the  pressure  was  130  kbars  and 
1600°C  for  experiments  at  300  kbars.  It  was  assumed  that  the  indicated 
temperature  closely  approximated  the  actual  temperature  rise  in  the  solder. 

This  interpretation,  which  assumes  no  fundamental  anomaly  of  the  Seebeck 
effect  in  the  shock  environment,  may  be  valid;  however,  the  ratio  of  the 
"measured  temperatures"  in  the  two  sets  of  experiments  is  too  small  compared 

with  the  ratio  of  pressures  for  any  reasonable  equation  of  state  for  solder. 

14 

Doran  and  Ahrens  reported  data  for  a series  of  experiments  in  which 

2 2 

the  junction  area  was  varied  from  0.2  mm  to  285  mm  and  the  materials  making 
up  the  junction  were  varied.  Signal  amplitudes  for  different  couples  were 
reported  to  be  in  the  same  ratio  as  their  zero  pressure  thermopowers,  but  a 
close  correlation  is  not  obvious  from  the  reported  data.  No  definite  relation- 
ship between  junction  area  and  amplitude  was  found.  A partially  successful 
diffusion-welding  process  of  the  junction  resulted  in  a lower  emf  in  one  of 
two  experiments.  The  authors  concluded  that  the  emf  is  probably  of  the  same 
origin  as  the  normal  thermoelectric  effect  and  that  it  is  more  nearly  propor- 
tional to  pressure  than  temperature.  In  these  experiments,  as  in  the  previous 
study  by  Jacquesson,  only  the  peak  amplitude  was  reported  and  compared  with 
the  expected  emf.  This  is  due  primarily  to  the  use  of  explosive  shock  gener- 
ation in  which  only  the  initial  pressure  is  well  known. 

1 5 

Palmer  and  Turner  conducted  a series  of  experiments  in  which  a 
constantan  wire  was  soldered  to  the  back  surface  of  a copper  plate  being 
impacted  from  the  front.  The  pressure  behind  the  shock  was  varied  from  40 
kbars  to  300  kbars.  The  voltage-time  profile  reported  shows  a 2 ysec  rise 
time  of  the  initial  voltage  and  a 20  ysec  pulse  with  random  excursions  of  up 
to  50%  of  the  peak  value.  The  authors  concluded  that  the  emf  was  a function 
of  final  shock  pressure  and  thereby  could  be  used  to  indicate  shock  strength 
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qualitatively.  It  is  not  obvious  that  the  reported  voltage  profile  can  be 
interpreted  in  this  way.  The  high  shock  pressure  in  a soldered  junction  of 
this  type  should  last  only  a few  hundredths  to  a few  tenths  of  a microsecond, 
depending  on  dimensions  and  configuration  of  the  solder  connection. 

Crosnier  ^al^.^^  reported  data  for  several  different  couples.  The 
junctions  were  formed  by  placing  a disc  of  one  material  on  the  back  surface 
of  the  other  with  no  radial  support,  or  by  placing  a rod  of  one  material  in 
an  insulated  hole  in  the  other.  The  authors  concluded  that  the  emf  was  due 
to  a voltage  source  of  low  internal  resistance,  and  that  the  junction  behaved 
qualitatively  like  a normal  thermocouple.  It  was  concluded  that  the  emf  was 
proportional  to  the  pressure  up  to  1600  kbars,  and  that  the  emf  was  of  such 
large  magnitude  that  no  classical  interpretation  was  possible.  It  was  pro- 
posed that  the  emf  could  be  due  to  the  formation  of  an  electronic  hot  gas 
with  a much  higher  temperature  than  the  crystal  lattice.  Although  these 
experiments  were  done  with  and  without  radial  pressure  support  behind  the 
shock  wave  in  the  active  element,  no  mention  is  made  of  the  effect  this  has 
on  the  observed  emf. 

In  1967,  Migault  and  Jacquesson^^  proposed  a theoretical  explanation 
for  the  abnormally  large  emfs  observed  in  their  previous  experimental  program. 
The  theory  is  based  on  the  radiation  pressure  exerted  on  the  conduction  elec- 
trons by  phonons  created  in  the  shock  front.  An  estimate  is  made  of  this 
effect  assuming  all  of  the  increase  in  internal  energy  is  in  the  form  of 
longitudinal  phonons  polarized  in  the  shock  propagation  direction.  This 
theory  is  said  to  explain  the  experimental  results  in  the  pressure  range  from 

80  kbars  to  300  kbars. 

1 8 

Conze  et  al^. , proposed  another  theoretical  explanation  of  the 
abnormally  high  emfs  in  this  type  of  experiment.  The  theory  is  based  on  an 

i 

t 

1 
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elevated  temperature  of  the  electrons  above  the  lattice  temperature  coupled 
with  an  electronic  imbalance  due  to  phonons  created  in  the  shock  front. 

Several  experiments  were  also  conducted,  including  one  in  which  mercury  was 
used  to  make  the  contact  in  the  junction.  Results  of  this  experiment  were 
reported  to  be  identical  to  those  in  which  the  surfaces  were  simply  pressed 

together.  The  one  reported  data  point  is  given  in  Table  2.1. 

19 

Buzhinskii  and  Samylov  conducted  a series  of  experiments  on  a copper- 
nickel  junction  made  from  disks  by  polishing  the  mating  surfaces  or  diffusion- 
welding the  disks  together.  No  distinction  was  made  in  the  reported  data 
between  these  two  types  of  junctions.  Results  are  given  in  the  table  for 
thirteen  separate  experiments  at  the  same  pressure.  The  diameter  of  the 
junction  was  varied  from  20  mm  to  32  mm  with  no  significant  change  in  the 
output.  The  mean  measured  temperature  was  reported  to  be  60%  greater  than 
the  calculated  temperature.  Although  a lot  of  scatter  is  present  in  the 
data,  which  may  be  partially  due  to  differences  in  shock  pressure,  the 
lowest  reported  emf  values  are  within  15%  of  the  expected  value.  No  indi- 
cation is  given  as  to  which  experiments  were  done  with  the  welded  junctions. 
These  results  suggest  that  the  higher  emf  values  observed  previously  may  be 
due  to  an  experimental  artifact  rather  than  to  a fundamental  phenomenon  such 

as  a high  electron  gas  temperature. 

20 

Lascar  and  Duage  conducted  a series  of  experiments  on  the  copper- 
constantan  junction  to  localize  the  source  of  the  abnormal  emf.  In  one 
series  of  experiments,  the  constantan  was  preheated  to  reduce  the  final  tem- 
perature difference  between  the  materials.  This  had  only  a "slight  effect" 
from  0°C  to  150°C  preheating.  In  another  series,  the  junction  was  diffusion 
welded  with  an  interdiffusion  layer  of  100  p and  an  alumina  ceramic  insulator 
was  used  around  the  copper  element,  which  was  inserted  in  a hole  in  the 
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constantan.  No  actual  data  are  given  for  these  experiments.  The  authors  do  | 

4 

not  draw  any  specific  conclusions  from  this  experiment,  but  the  qualitative  | 

sketch  of  the  voltage-time  profile  shows  a relatively  small  initial  step  in  | 

emf  followed  by  a ramping.  Without  more  specific  data,  an  interpretation  of  | 

the  significance  of  these  results  in  terms  of  our  observations  is  impossible.  I 

i 

The  authors  conclude  that  the  abnormal  emf  is  not  localized  at  the  junction  ] 

or  in  the  shock  front  and  therefore,  the  emf  source  must  be  distributed  in  | 

the  bulk  of  the  materials  making  up  the  junction.  This  conclusion  does  not  \ 

appear  to  be  justified  by  the  qualitative  results  reported. 

21  ' 
In  1970  Jacquesson  ^ , reviewed  the  results  of  previous  experi-  | 

ments  in  this  area,  and  reiterated  the  conclusion  that  the  emf  has  a linear  3 

I 

pressure  relation  in  this  type  of  experiment  in  the  form  E(mv)  = 0.145P  ! 

i 

(kbars)  +11.  No  reference  was  made  to  the  generally  lower  results  of  j 

Buzhinskii  and  Samylovr  Some  new  experiments  were  conducted  in  which  the  i 

plane  of  the  advancing  shock  front  was  perpendicular  to  the  junction  plane.  ’ 

In  these  experiments,  the  emf  increased  monotonically  with  time  and  reached 
a maximum  when  the  entire  junction  was  pressurized.  The  final  emf  was 
reported  to  be  the  same  as  in  the  case  where  the  shock  front  and  junction 
plane  are  parallel . 

22 

The  results  by  Mineev  et  al_. . on  the  lanthanides  ytterbium,  europium 
and  cerium  are  included  in  Table  2.1.  The  junctions  were  made  by  pressing 
1-2  cm  disks  of  the  sample  on  the  back  of  an  aluminum  buffer  plate.  The 
authors  conclude  that  the  measured  emf  is  not  connected  with  contact  poten- 
tial effects  and  is  principally  due  to  volume  redistribution  of  charge  in 
the  shock  front. 

23 

Bordzilovskii  et  ^. , conducted  a series  of  experiments  in  which 


the  angle  between  the  plane  of  the  junction  and  the  plane  of  the  shock  wave 


was  varied.  The  rise  time  of  the  measured  emf  was  equal  to  the  time  neces- 
sary for  the  shock  wave  to  sweep  across  the  junction.  The  authors  conclude 
that  the  emf  is  due  to  an  internal  contact  potential  difference  mechanism. 

The  data  reported  for  copper-nickel  and  duralumin-nickel,  where  the  incidence 

24 

angle  was  zero,  are  included  in  Table  2.1.  In  another  set  of  experiments 

using  an  external  shunt  of  low  impedance,  the  authors  concluded  that  internal 

-2  2 
resistance  of  this  emf  source  was  about  10  for  a contact  area  of  78.5  mm 

and  was  inversely  proportional  to  the  area.  The  objective  of  these  studies 

was  to  characterize  the  electrical  response  of  a junction  in  order  to  use  it 

as  a pressure-time  profile  transducer.  No  attempts  were  made  to  determine 

the  cause  of  the  large  emf  as  compared  with  the  normal  thermoelectric 

response. 

25 

In  1972  Migault  et  aj[. , gave  a modified  theoretical  treatment  of 

the  response  of  a bimetallic  junction  to  shock  compression.  The  theory  takes 

account  of  the  radiation  pressure  of  the  phonon  gas  generated  by  the  shock. 

The  resultant  effect  is  found  by  a resolution  of  a coupled  set  of  linearized 

Boltzmann  equations  for  phonon  and  electron  distributions.  The  proposed 

theory  purportedly  explains  the  high  observed  emfs  in  the  pressure  range 

from  0 to  300  kbars.  Again,  no  reference  is  made  to  the  lower  emf  values 

observed  by  Buzhinskii  and  Samylov,  which  would  not  agree  with  this  theory. 

26 

In  1974,  Bordzilovskii  ^ , reported  using  a copper-constantan 

junction  to  measure  a pressure  profile  in  copper.  The  pressure  wave  in 
copper  was  produced  by  impacting  a flyer  plate  on  an  intermediate  layer  con- 
sisting of  polystyrene  foam,  liquid  nitrogen  or  liquid  hydrogen.  The  junc- 
tion was  made  by  pressing  a constantan  disk  of  6 mm  diameter  on  the  back 
surface  of  the  copper  plate.  The  emf-pressure  calibration  curve  which  was 
used  is  given  in  Table  2.1  and  was  apparently  measured  in  a separate  set  of 


{ 


experiments.  The  authors  conclude  that  the  recording  method  is  satisfactory 
to  within  10-15%.  However,  the  authors  did  note  a discrepancy  of  as  much  as 
65%  between  the  calculated  and  measured  final  pressures.  The  reported 
voltage-time  profiles  do  not  appear  to  be  consistent  with  the  conclusion  that 

a direct  measurement  of  pressure  was  being  made. 

27  2R 

Nesterenko  e;t  al^- . and  Nesterenko  conducted  a series  of  experi- 
ments using  various  pairs  of  metals.  In  one  series  on  copper-nickel  at 
400  kbars,  half  of  the  experiments  were  done  with  emery-polished  surfaces  and 
I the  other  half  were  done  using  "rough"  surfaces.  The  polished  surfaces 

resulted  in  40%  reduction  in  amplitude.  The  authors  conclude  that  the  high 
emf  measured  in  this  type  of  experiment  is  due  to  a high  nonequilibrium  con- 
I tact  temperature. 

In  summary,  the  results  of  previous  studies  in  this  area  are  con- 
fusing, but  some  conclusions  can  be  drawn.  The  emf  is  of  the  appropriate 
sign  to  be  interpreted  on  the  basis  of  a normal  thermoelectric  response  to  an 
elevated  temperature.  The  emf  observed  in  every  case  is  higher  than  predicted 

on  a thermoelectric  basis,  but  the  discrepancy  varies  from  as  little  as  15% 

r 19  28 

in  some  experiments  by  Buzhinskii  ^ , and  Nesterenko  to  as  much  as  a 

factor  of  twenty  or  more.  This  also  indicates  the  general  nonreproducibility 

which  is  evident  both  among  different  investigations  and  within  any  particular 

investigation,  Nesterenko  has  shown  that  when  polished  surfaces  are  used  at 

the  interface,  the  resultant  emf  is  lower  than  when  rough  surfaces  are  used. 

The  results  of  diffusion-welding  the  junction  are  not  conclusive  but  suggest 

a lower  resultant  emf.  Signals  are  generally  not  steady  in  time  and  are 

usually  noisy.  The  investigations  which  have  been  undertaken  up  to  this 

point  have  been  directed  toward  either  the  characterization  of  this  effect 

to  make  it  useful  as  a pressure  transducer  or  toward  the  resolution  of  the 
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conflict  between  the  predicted  and  observed  results.  In  the  first  case  the 
nonreproducibility  among  investigators  puts  the  results  in  doubt;  and  in  the 
second  case,  the  results  and  conclusions  among  authors  vary  widely  and  no 
general  resolution  of  the  problem  exists. 
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CHAPTER  3 


SOURCES  OF  ANOMALOUS  EMF 

In  the  first  phase  of  the  present  program  an  experimental  configu- 
ration was  developed  in  which  the  measured  emf  could  be  predicted  by  direct 
comparison  with  the  idealized  experiment  of  Chapter  1.  This  was  done  by 
recognizing  and  eliminating  any  mechanism  that  would  cause  an  emf  to  be  mea- 
sured which  was  not  associated  with  the  junction  response  to  shock  compres- 
sion as  discussed  in  Chapter  1.  Only  the  fruits  of  this  effort,  specifically 
the  three  mechanisms  which  were  isolated,  will  be  described  in  this  chapter. 
The  details  of  the  associated  experiments  will  be  included  in  Appendices  C, 

D,  and  E.  The  final  configuration  will  be  described  in  Chapter  4. 

3.1  Electrical  Noise  Observed  When  Circuit  Is  Completed  during 
Shock  Experiment. 

We  have  observed  that  when  a low  impedance  circuit  (~0.1  ohm)  is 
completed  during  a shock  experiment  electrical  noise  is  observed  in  the  cir- 
cuit. In  a thermoelectric  experiment  the  two  possibilities  for  this  to  be 
observed  are  when  part  of  the  eventual  circuit  is  in  the  projectile,  or  when 
a small  gap  between  components  is  closed  by  a shock  wave.  In  either  case  the 

resulting  noise  interferes  with  the  measurement  of  the  thermoelectric  emf. 

19 

Buzhinskii  and  Samylov  observed  that  when  a shock  wave  passes 
through  a series  of  disks  that  are  nominally  in  contact,  and  electrically  in 
series  with  the  recording  instrument,  a noise  spike  is  observed  for  each 
interface.  The  details  of  the  experiment  and  the  amplitude  of  the  noise 
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were  not  given.  In  much  of  the  previous  work  (see  Chapter  2)  the  thermo- 
electric circuit  was  prepared  by  simply  holding  the  elements  in  contact.  In 
studies  where  only  the  initial  peak  signals  were  reported,  the  values  may  be 
anomalous  because  of  the  simultaneous  noise. 

In  experiments  done  here,  this  type  of  noise  was  observed  when 
attempts  were  made  to  complete  the  circuit  at  impact.  The  objective  was  to 
have  one  part  of  the  bimetallic  pair  in  the  projectile  and  thus  produce  sym- 
metric shock  waves  in  the  two  metals  as  in  the  idealized  experiment  of 
Chapter  1.  Several  experiments  were  conducted  to  investigate  the  problem  and 
the  details  are  given  in  Appendix  C.  Although  the  fundamental  source  of  this 
noise  is  not  well  understood,  we  consistently  observed  a non-zero  signal  from 
as  much  as  5 psec  before  impact.  At  impact  a positive  and/or  negative  spike 
greater  than  100  mv  was  observed  and  a high  frequency  (>50  mHz)  noise  per- 
sisted for  up  to  0.5  psec  after  impact.  Initially  our  efforts  were  directed 
toward  minimizing  or  eliminating  this  problem  using  electrical  and  magnetic 
shielding  techniques,  but  we  eventually  concluded  that  closure  of  the  circuit 
at  impact  must  be  avoided  completely  to  eliminate  this  problem. 

3.2  Effect  of  Shock  Demagnetization. 

In  several  earlier  studies  involving  the  shock  loading  of  bimetallic 
junctions,  a ferromagnetic  material,  iron  or  nickel,  was  used  as  one  of  the 
active  elements. It  has  been  shown  in  numerous  investi- 
gations that  ferromagnetic  materials  undergo  a demagnetization  process  upon 
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the  passage  of  a strong  shock  wave.  It  is  necessary,  therefore,  to 

determine  what  effect  the  change  in  magnetic  state  has  on  the  results  of  such 
experiments. 

The  following  mechanisms  for  shock  induced  demagnetization  have  been 
identified:^^ 


1.  First-order  phase  transitions  from  ferromagnetic  to 
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nonferromagnetic  phases  such  as  in  iron. 

2.  Second-order  phase  transitions  between  ordered  and  dis- 
ordered states  in  combination  with  a lowering  of  the 
Curie  temperature  and  shock  heating. 

3.  Strain  induced  magnetic  anisotropy  in  which  the  direction 
of  shock  propagation  becomes  an  easy  axis.  This  leads  to 
an  effective  demagnetization  in  the  transverse  direction. 

Iron  was  used  in  only  two  previous  investigations^^’^®  and  was  not 
given  primary  consideration  in  either  study.  Since  no  work  was  done  here 
involving  iron,  attention  will  be  focused  on  the  shock  demagnetization  of 
nickel.  There  are  no  known  first-order  phase  transitions  in  nickel  below 
500  kbars  thus  eliminating  the  first  possible  mechanism. 

Experimental  evidence  of  some  form  of  demagnetization  process  in 
34 

nickel  was  seen  by  Wong  at  23  kbars  in  the  form  of  an  eddy  current  spike 
in  a resistivity  measurement.  It  would  not  be  possible  to  conclude  from  this 
type  of  experiment  which  of  the  two  mechanisms  listed  above  is  active. 

Although  no  study  of  the  third  mechanism  has  been  done  using  nickel, 

a realignment  of  the  magnetization  along  the  shock  propagation  direction  has 

been  observed  for  other  materials,  including  nickel  ferrite  and  yttrium  iron 
29,30 

garnet. 


Grady  ei^  aj^.  have  shown  for  polycrystalline  materials  that  the 
direction  of  magnetization  is  determined  by  a minimizing  of  the  internal 
energy  expressed  as 


E = -H„*M  + Besin^e 
e 


(3.1) 


where  E is  the  internal  energy,  is  the  external  magnetic  field,  M is 


r—  ^ 

44 

the  magnetization,  e is  the  uniaxial  strain,  0 is  the  angle  between  the 
direction  of  the  applied  field  and  the  direction  of  the  magnetization,  and 
B is  a constant  defined  as 

B = (-3/5)[(Cii  - Ci2)Yiqq  3C^^Yiii]  ' (3.2) 

Here  C..  and  y. ..  are  the  elastic  and  magnetostrictive  constants,  respec- 

1 J 1 J K 

tively.  The  internal  energy,  E,  is  a sum  of  the  mechanical  elastic  energy, 
the  magnetoelastic  energy,  and  the  energy  due  to  interaction  of  the  external 

■ 

magnetic  field  and  the  magnetization.  ] 


The  experiments  to  measure  this  effect  are  done  by  magnetizing  the 
sample  in  a direction  perpendicular  to  the  shock  propagation  direction  with  I 

i 

: an  external  field.  When  the  shock  propagates  through  the  material  the  magne- 

j tization  direction  changes  to  minimize  the  energy  in  Eq.  (3.1).  When  B is  i 

t i 

positive,  as  is  the  case  for  nickel  ferrite,  YIG,  and  nickel,  the  second  term  i 

in  Eq.  (3.1)  is  minimized  when  0 is!90°,  thus  tending  to  realign  the  magne- 

tijzation  direction  with  the  shock  propagation  direction.  This  is  observed 
I as  a demagnetization  in  the  transverse  direction.  In  the  special  case  where 

i 

I no  external  magnetic  field  is  applied,  a complete  realignment  is  expected 

I along  the  shock  propagation  direction. 

I A series  of  8 experiments  was  conducted  in  our  study  of  the  demagne- 

I tization  process  and  its  effect  on  bimetallic  junction  experiments.*  The 

details  of  these  experiments  are  given  in  Appendix  D.  The  experimental  pro- 
gram was  initially  undertaken  to  determine  whether  or  not  shock  compression 

*The  term  "bimetallic  junction  experiment"  is  in  reference  to  the 
type  of  experiment  that  is  indicated  schematically  as  the  idealized  experi- 
ment of  Chapter  1.  This  could  include  a number  of  specific  configurations, 
but  each  would  have  in  common  a closed  circuit  of  low  internal  impedance 
(-O.Ol  fi),  two  metals  to  form  the  junction  with  a shock  wave  traveling  in 
each,  and  a recording  system  sensitive  to  the  small  expected  emfs 

(0.1  - 50  mv).  ! 
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of  a bimetallic  junction  in  which  both  metals  are  the  same  would  produce  a 
null  result.  A null  result  is  predicted  by  every  previously  proposed  theory. 

In  the  first  such  experiment  (76-032)  the  recorded  emf  was  not  zero.  At  this 
point  we  proposed  several  possible  causes  for  the  non-zero  signals,  and  the 
next  three  shots  (76-037,  76-039,  76-040)  were  conducted  before  we  realized 
that  the  emf  was  probably  due  to  demagnetization. 

Three  shots  76-051,  76-054,*  and  76-062  were  done  specifically  to 
determine  whether  or  not  a demagnetization  of  the  nickel  did  occur  under  the 
experimental  conditions  of  interest,  and  if  this  process  would  interact  with 
the  recording  circuit  and  result  in  the  characteristic  signals.  In  shots 
76-051  and  76-054  the  magnetization  of  the  nickel  was  controlled.  In  76-051  ' 
the  nickel  was  used  as  received  and  had  a magnetic  field  just  outside  the 
body  of  about  0.8  Gauss  in  directions  parallel  and  perpendicular  to  the  shock 
propagation  direction.  In  76-054  the  nickel  was  demagnetized  with  a field 
just  outside  the  body  of  less  than  0.05  Gauss  in  any  direction.  The  recorded 
emf  was  about  -40  mv  in  the  first  of  the  shots  and  about  -7  mv  in  the  second. 
(See  Fig.  D.5.)  Therefore  a reduction  of  the  pre-shot  magnetization  in  the 
nickel  by  a factor  of  16  reduced  the  signal  amplitude  by  about  a factor  of  6. 
This  effect  of  the  initial  magnetization  would  be  expected  if  the  signals 
were  caused  by  a shock  demagnetization. 

In  another  shot,  76-062,  the  active  element  was  changed  from  nickel  to 

j. 

Ni chrome  V'  (80%  nickel  2® chromium)  which  is  not  ferromagnetic.  The  experiment 
was  the  same  in  all  other  respects  as  shot  76-032.  In  shot  76-032  an  emf  of 
less  than  -100 mv  was  observed  and  in  76-062  no  emf  was  observed  during  the  time 

*The  author  wishes  to  acknowledge  the  efforts  of  Dr.  Jerry  Dick  who 
guided  the  work  in  this  area  and  conducted  shots  76-051  and  76-054. 

Trademark  of  the  Driver  Harris,  Co.,  Harrison,  New  Jersey. 


of  shock  transit  through  the  target  to  within  the  recording  capability  of 
the  instrumentation  system  of  0.01  mv.  (See  Fig.  D.l  and  D.5.) 


46 


Results  of  these  experiments  show  conclusively  that  an  emf  is  obser- 
ved in  this  type  of  experiment  whefi  one  of  the  active  elements  is  ferro- 
magnetic. The  emf  is  due  to  shock  demagnetization  of  that  element.  This 
puts  the  results  of  any  previous  work  done  with  nickel  or  iron  in  doubt.  It 
also  precludes  the  use  of  ferromagnetic  materials  in  this  type  of  experiment 
unless  this  effect  can  be  completely  nullified  in  some  other  way. 

3.3  Effect  of  Surface  Layer  between  Elements. 

In  most  of  the  work  previously  reported  in  which  the  response  of  a 
bimetallic  junction  to  shock  loading  has  been  investigated,  the  observed  emf 
is  2 to  15  times  that  predicted  assuming  a thermoelectric  source  and  reason- 
able shock  compression  temperatures.  One  possible  contributing  factor  to 
these  erratic  but  consistently  high  results  is  that  the  interface  between  the 
two  materials  reaches  a much  higher  temperature  than  the  bulk  material.  If 
indeed  a thermoelectric  source  is  responsible  for  the  emf  recorded,  the 
inferred  temperature  would  be  the  local  temperature  in  the  immediate  vicinity 
of  the  boundary  between  the  two  materials.  In  this  case,  the  temperature 
rise  of  the  interface  above  the  bulk  shock  compression  temperature  would  be 

zero  only  if  a perfect  plane  separated  the  two  materials.  Urtiew  and 
35 

Grover  have  studied  the  possible  distribution  of  temperature  in  both  dis- 
tance and  time  near  an  interface  between  two  materials  under  shock  loading 
conditions.  Two  models  were  proposed,  one  in  which  the  intermediate  zone 
was  characterized  by  a gap  as  might  be  the  case  for  highly  polished  surfaces 
and  another  in  which  the  intermediate  zone  is  considered  as  a porous  solid. 
Evidence  is  given  for  the  acceptance  of  the  porous  material  model  as  the  more 
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correct  even  for  highly  polished  surfaces  for  which  the  interface  zone  is  on 
the  order  of  one  micron  thick. 

It  was  shown  in  experiments  involving  heat  generation  due  to  the 

36 

friction  of  sliding  by  Bowden  and  Tabor  that  hot  spots  on  the  order  of  a 
micron  or  less  would  produce  an  emf  in  a thermoelectric  circuit.  A diagram 
of  the  type  of  experiments  conducted  is  seen  in  Fig.  3.1a.  In  Fig.  3.1b, 
the  surfaces  in  contact  are  shown  microscopically.  It  was  proposed  that  in 
the  small  regions  which  make  contact,  a very  high  temperature  was  produced 
due  to  friction,  and  the  purpose  of  the  experiments  was  to  determine  if  this 
was  true  by  measuring  the  thermoelectric  "response  of  a constantan  rod  rubbing 
on  a copper  disk.  Transient  temperatures  on  the  order  of  1000°C  were  mea- 
sured for  time  of  less  than  one  millisecond.  These  results  indicate  that 
the  thermoelectric  circuit  involved  in  the  shock  experiments  would  respond 
to  a region  of  high  temperature  on  the  order  of  several  microns  at  the  junc- 
tion interface. 

Consider  a plane  interface  between  copper  and  constantan  with  the 
following  properties: 

1.  The  intermediate  region  will  be  treated  as  a porous 
layer  of  thickness  2a. 

2.  The  porous  material  will  have  the  properties  of  copper 
foam. 

3.  Complete  compaction  is  achieved  in  the  first  shock 
transit. 

Figures  3.2  and  3.3  show  that  the  pressure-particle  velocity  state  reached 
in  the  porous  layer  after  one  reflection  may  be  nearly  the  same  as  the 
initial  state  in  the  solid.  However,  the  temperature  in  the  interface 
region  is  higher  than  that  of  the  bulk  material. 


Constantan 


(a) 


Fig.  3, 1 . --inustration  of  Experiments  by  Bowden  and  Tabor  on 
Friction,  (n  (a)  an  experiment  to  measure  the  temper- 
ature produced  locally  by  friction  is  illustrated,  and 
in  (b)  a microscopic  view  of  the  sliding  surfaces  is 
shown. 
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Up 


Fig.  3. 2. --Pressure-Particle  Velocity  Diagram  for  Shock  Wave 


Interactions  in  a Porous  Layer  between  Solid 
Regions. 


The  temperature  in  the  porous  material  can  be  estimated  using  a 
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Christian  and  Walsh  technique  (see  section  1.3),  in  which  the  Hugoniot  of 
the  porous  material  is  characterized  before  compaction  by  the  relation 

P = A - bx  (3.3) 

where  P is  the  pressure  and  x is  the  compression  defined  by 

_ present  volume 

^ ~ full  density  volume  at  zero  pressure 


Here  A and  b are  positive  constants.  This  straight  line  pressure-volume 
curve  is  used  until  it  intersects  the  principal  Hugoniot  of  the  full  density 
material,  which  is  assumed  to  occur  when  compaction  to  full  density  is  com- 
plete. A more  nearly  exact  relation  for  P(v)  can  be  obtained  using  the 
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P-a  equation  of  state  of  Herrman.  However,  for  temperature  calculations 
the  difference  is  less  than  10%  up  to  500  kbars  and  the  P-a  form  is  much  more 
difficult  to  manipulate.  It  is  also  assumed  that  C^,  y,  and  p are  the 
same  as  for  the  zero  pressure  full  density  material  at  the  point  where  com- 
paction is  complete. 

After  compaction  is  complete,  the  Hugoniot  for  the  porous  layer  can 

be  calculated  from  the  full  density  Hugoniot  by  using  the  Griineisen  equation 

of  state  to  take  into  account  the  added  internal  energy  due  to  irreversible 
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work  done  in  compaction  (see  Boade).  Consider  two  independent  shock  situ- 
ations for  the  porous  and  full  density  materials.  To  each  of  these  single 
shocks  one  may  apply  the  Rankine-Hugoniot  equation 


E 


S 


^Pp(v 


^Ps(v 


Po 

So 


(3.4) 


(3.5) 


where  the  subscripts  P and  S refer  to  porous  and  full  density  material, 
respectively.  Now  to  find  Pp(v)  we  use  the  Gruneisen  equation 

Pp  - Pj  = (Y/v)(Ep  - E3)  (3.6) 

which  is  just  a relation  between  two  equilibrium  states  of  the  same  material 
at  constant  volume.  Substituting  from  Eqs.  (3.4)  and  (3.5),  Eq.  (3.6) 
becomes 


^P  “ 2 ^^P^Vo  “ ''P^  ■ ^S^''So  " 


where  it  is  assumed  Ep^  = E^^.  To  find  the  relation  for  Pp(v)  for  a 
given  which  is  known,  let  Vp  = v^  = v,  so  that 


P 


P 


(3.8) 


where  y/v  = G = constant.  With  the  Hugoniot  just  constructed,  it  is  pos- 
sible to  estimate  the  temperature  in  the  porous  layer  by  a numerical  inte- 
gration method  outlined  in  section  1.3.  Temperature  rise  in  a single  shock 
to  various  pressures  and  from  various  initial  densities  is  given  in  Table  3.1 
for  copper.  These  temperatures  are  as  much  as  a factor  of  20  higher  than 
shock  compression  temperatures  in  the  solid.  It  should  also  be  noted  that 
the  temperatures  in  the  porous  layer  are  roughly  proportional  to  the  shock 
pressure  for  a given  initial  porosity. 

The  temperature  profile  just  after  the  shock  wave  passes  through  the 
interface  is  shown  in  Fig.  3.4.  Evolution  of  this  initial  profile  in  time 
must  be  considered  in  terms  of  the  relative  time  for  decay  of  the  temperature 
profile  and  the  lifetime  of  the  experiment.  The  maxinum  temperature,  which 
will  always  be  at  y = 0,  decreases  with  time.  Depending  on  the  relative  time 


dependence,  one  would  expect  either  to  see  a predictable  rate  of  decay  of 
the  thermal  emf,  or  in  the  case  of  long  decay  times,  the  only  recorded 
thermal  emf  would  be  that  due  to  the  initial  temperature  profile. 
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TABLE  3.1 .--Temperature  Rise  in  Shocked  Porous  Copper 


Pressure 

(kbar) 

1.5 

Temperature  Rise 
Initial  Volume  v/v^* 
1.3  1.2  1.1 

1.0** 

150 

1500 

920 

635 

345 

76 

200 

1925 

1200 

830 

465 

109 

250 

2360 

1480 

1025 

585 

149 

300 

2800 

1770 

1235 

715 

196 

350 

3240 

2065 

1460 

850 

249 

400 

3715 

2375 

1685 

990 

309 

450 

4150 

2675 

1915 

1140 

376 

500 

4630 

3000 

2155 

1300 

449 

* 

''s 

is  the 

specific  volume  for 

solid 

copper. 

**Values  of  T for  full  density  copper. 

Values  used  for  the  constants  are: 

Y = 1.99 

= 8.93  gram/cm^  (full  density) 

= 3.718  X 106  ergs/g°K 
Cq  = 3.94  mm/ysec 
S = 1.489 

where  + S Up  for  full  density  material. 

The  time  dependence  of  the  initial  distribution  can  be  studied  using 


the  one  dimensional  heat  flow  equation.  It  can  be  shown  (see  Appendix  E) 
that  the  temperature  at  y = 0,  T(0,t),  is 


T takes  on  the  values  0.01  usee,  0.25  usee,  and  1 usee  for  porous  layers  of 
1 U,  5 u,  and  10  u,  respeetively . A plot  of  T(0,t)/T2  versus  t is  given 
in  Fig.  3.5.  Typieal  experiment  life  times  are  about  1 usee  in  most  previous 
work.  Therefore,  if  an  effeetive  porous  layer  of  10  u existed  at  the  inter- 
face, the  temperature  would  decay  less  than  40%  during  the  experiment. 

To  show  the  importance  of  the  temperature  at  y = 0,  consider  an 
experiment  in  which  plane  impact  occurs  between  slabs  of  copper  and  constantan. 
The  receding  shock  in  constantan  and  the  advancing  shock  in  copper  are  indi- 
cated in  Fig.  3.6  along  with  a temperature  profile  which  includes  a high 
temperature  region  about  the  origin.  Using  Eq.  (1.2b)  the  emf  measured 
between  a and  b in  this  circuit  is  given  by 

^2  Tq 

<^a-b  = - I ^Cu  - j dT  . (3.10) 

^0  "'^2 

If  we  assume  for  simplicity  that  and  are  independent  of  temper- 

ature and  pressure  this  becomes 

‘*’a-b  ^ “^Cu^^2  “ ^0^  ■ ^Co^^O  ‘ ^2^ 


“^a-b  ""  ^^Co  ■ ^Cu^^^2  ■ ^0^  ‘ (3.11) 

The  emf  observed  in  such  a circuit  is  independent  of  T.|  and  T^,  which 
represent  the  shock  compression  temperatures,  and  is  a function  only  of  T^, 
which  is  T(0,t)  given  in  Eq.  (3.9).  This  is  true  even  if  T^  is  only 
slightly  higher  than  the  shock  compression  temperature. 


Distance 


Fig.  3.6. — Thermoelectric  Circuit  and  Associated  Temperature 
Profile  with  Porous  Interface  Layer. 


Two  experiments  were  conducted  to  investigate  the  effect  of  surface 
preparation  on  measured  emf  in  a copper-constantan  junction  experiment.  The 
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experiments  are  discussed  in  detail  in  Appendix  E.  In  one  experiment  the 
interface  was  prepared  by  lapping  the  two  surfaces  and  in  the  other  the  sur- 
faces were  highly  polished.  Measured  emfs  for  the  lapped  and  polished  sur- 
face shots  were  -19  ± 1 mv  and  -16  ± 2 mv,  respedtively.  The  voltage-time 

I 

profiles  are  given  in  Fig.  E.2.  The  predicted  efhf  based  on  a shock  compres- 

\ 

sion  temperature  of  the  bulk  materials  is  about '-3  mv.  This  indicates  in 
terms  of  this  model  that  both  the  lapped  and  polished  surfaces  result  in  a 

i 

porous  region,  but  that  polished  surfaces  result  in  a lower  effective  porosity. 
Since  only  two  shots  were  done  no  firm  conclusions  can  be  drawn. 

In  many  of  the  experiments  previously  reported  in  which  an  anomalously 
high  emf  was  recorded,  the  surfaces  involved  are  not  well  characterized.  A 
complete  understanding,  therefore,  of  the  relationship  of  these  results  to 
the  problem  of  surface  heating  is  not  possible.  In  experiments  conducted 
here,  it  was  noted  that  even  for  highly  polished  surfaces,  a wave-like  dis- 
turbance develops  on  the  polished  surface  with  a thickness  on  the  order  of  a 
micron.  This  indicates  that  though  careful  surface  preparation  may  reduce 
the  problem,  it  would  not  be  possible  to  completely  eliminate  these  effects. 
Bowden  and  Tabor  have  shown  through  surface  resistance  measurements  that 
the  actual  area  in  contact  between  two  surfaces  can  be  more  than  an  order  of 
magnitude  less  than  the  geometrical  area  even  for  highly  polished  surfaces. 

This  too  indicates  that  even  under  the  best  of  conditions,  a region  in  which 
there  are  voids  must  exist  between  the  two  bulk  materials  of  interest. 


Results  of  our  experiments  indicate  that  polished  surfaces  at  the 
junction  may  result  in  a lower  emf  than  when  rough  surfaces  are  used.  We 
have  shown  that  the  "measured  temperature,"  or  the  temperature  to  which  the 
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[ junction  as  a thermocouple  responds,  is  the  temperature  in  the  immediate 

E 

^ vicinity  of  the  junction,  and  is  independent  of  the  temperature  behind  the 

shock  front.  Finally,  we  have  shown  that  if  the  interface  region  is  treated 
as  a porous  solid,  the  temperatures  reached  by  a shock  compression  are  high 
I . enough  to  explain  the  anomalously  large  emfs  observed  previously.  It  seems 

quite  reasonable  to  conclude  that  the  measured  response  to  shock  loading  of  a 
bimetallic  junction,  in  which  the  specimens  are  simply  held  in  contact,  is 
partially  or  perhaps  totally  due  to  the  conditions  at  the  junction  interface. 


i 


i 


CHAPTER  4 


EXPERIMENTAL  TECHNIQUE 

The  experimental  configuration  and  techniques  described  here  were 
developed  to  avoid  the  problems  discussed  in  Chapter  3 and  still  retain  the 
essential  conditions  of  the  idealized  experiment  of  Chapter  1.  The  actual 
experiments  will  be  described  first  followed  by  a discussion  of  how  the 
specific  problems  were  avoided.  Finally,  a comparison  will  be  made  between 
the  real  experiments  and  the  idealized  experiment. 

4.1  Experimental  Configuration. 

The  basic  impactor  and  target  assemblies  are  illustrated  in  Fig.  4.1. 
The  constantan  is  diffusion-welded  to  the  copper  buffer  as  shown,  with  a 
ceramic  and  copper  guard  ring  surrounding  it  for  lateral  support.  The  ceramic 
was  machined  to  a close  tolerance,  and  fitted  between  the  copper  and  con- 
stantan, by  a process  described  in  section  4.4.  The  copper  leads  are 
attached  as  indicated  with  solder.  The  copper  guard  ring  is  in  intimate  con- 
tact with  the  buffer  plate,  but  is  not  involved  in  the  thermoelectric  circuit. 
Nominal  dimensions  are  listed  in  Fig.  4.1  for  each  of  the  components. 

Shock  pressure  is  generated  by  plane  high  velocity  impact  of  a copper 
or  tungsten  alloy  impactor  on  the  copper  buffer.  The  plane  shock  wave  travels 
through  the  copper  buffer,  and  uniformly  loads  the  copper-constantan  junction. 
A shock  wave  also  travels  out  radially  from  the  edge  of  the  impactor  toward 
the  outer  edge  of  the  copper  buffer.  The  dimensions  here  were  chosen  so  that 
no  pressure  disturbance  reaches  the  point  on  the  buffer  where  the  measurement 
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lead  is  attached  before  measurement  of  the  thermoelectric  emf  is  completed. 

The  impactor  is  thick  enough  that  the  pressure  at  the  impact  plane  is  main- 
tained until  the  shock  wave  in  the  constantan  reaches  the  free  surface. 

There  will  be  a slight  perturbation  of  the  pressure  profile  in  the 
constantan  due  to  the  ceramic.  The  shock  impedance,  PqU^,  of  this  type  of 
ceramic,  which  will  determine  the  final  pressure  state  due  to  a single  shock, 
is  essentially  the  same  as  that  of  copper  below  the  Hugoniot  elastic  limit  of 
the  ceramic.  However,  in  the  plastic  region  above  80  kbars,  a more  compli- 
cated situation  arises  and  the  perturbation  might  become  quite  large  at  high 
pressures.  This  effect  will  be  discussed  further  when  the  results  of  the 
high  pressure  shots  are  presented  in  Chapter  5.  The  shock  velocity  in  the 
ceramic  is  greater  than  the  shock  velocity  in  constantan.  Therefore,  the 
wave  running  ahead  might  be  expected  to  cause  some  problems,  but  none  have 
been  identified. 

4.2  Materials. 

The  four  basic  materials  from  which  the  target  and  impactor  are  con- 
structed are  copper,  constantan,  ceramic  for  the  guard  ring,  and  a tungsten 
alloy  used  as  the  impactor  for  the  high  pressure  experiments.  The  copper 
used  throughout  was  oxygen-free,  high  conductivity  copper,  usually  designated 
OFHC  copper. 

The  constantan  used  in  these  experiments  was  obtained  from  two  sources 
The  material  designated  A-constantan  was  manufactured  in  a polycrystalline 
form  from  45  wt  % nickel  and  55  wt  % copper.  The  starting  materials  were 
99.99%  pure.  In  the  process  of  sectioning  the  material,  it  was  noted  that 
voids  of  approximately  one  micron  diameter  were  occasionally  visible.  No 
more  than  10-15  pores  of  this  type  were  seen  on  a single  section  of  one  inch 
diameter.  Densities  measured  for  this  material  were  within  the  measurement 
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error  of  the  expected  value  of  8.9  gram/ cm  . More  information  for  this  and 
the  following  materials  is  included  in  Table  4.1. 


TABLE  4.1  .—Materials  Used 


Supplier 

Material  Type 

Identification 
Used  Here 

Analyzed 

Composition 

Metals  Research* 

45  wt%  Nickel 

55  wt%  Copper 
Alloy 

A-Constantan 

44.4  wt%  to 

46.9  wt%  Nickel 

U.S.  Bureau  of 

Mi nes+ 

44.5  wt%  Nickel 

55.5  wt%  Copper 

1%  Manganese 

Alloy 

B-Constantan 

54.2%  Copper 
44.2%  Nickel 

1%  Manganese 
0.3%  Carbon 

Teledyne** 

17-D  Densalloy 

Tungsten  Alloy 

90.91%  Tungsten 
4%  Iron 

3%  Copper 

2%  Nickel 

♦Metals  Research  Ltd.  Mel  bourn,  Royston,  Herts  SG8  6EJ,  England. 


U.S.  Bureau  of  Mines,  Albany,  Oregon. 

**Teledyne  Powder  Alloys,  Clifton,  New  Jersey. 

The  material  designated  B-constantan  was  supplied  by  the  U.S.  Bureau  of 
Mines.  This  material  was  manufactured  in  a polycrystalline  form  from  54.5 
wt%  copper,  44.5  wt%  nickel,  and  1 wt%  manganese.  The  starting  materials 
were  commercially  pure,  but  the  alloying  was  done  in  a carbon  vessel  and 
later  analysis  shows  0.3%  carbon  exists  in  the  material.  The  density  of  this 

3 

material  measured  8.82  ± 0.05  g/cm  . 

The  ceramic  used  for  the  guard  ring  was  supplied  by  Coors  Porcelain 
Company,  and  is  designated  AD-998  Alumina.  The  ceramic  was  supplied  as  a 
tubing  and  is  nominally  99.8%  AI2O2. 

The  tungsten  alloy  which  was  used  as  the  impactor  for  the  higher 
pressure  shots  was  supplied  from  two  different  production  lots.  The 


properties  measured  such  as  density,  alloy  composition,  and  yield  strength, 
are  not  significantly  different  between  lots.  Hugoniot  data  for  this  material 
are  given  in  Appendix  B.  The  concentration  of  the  various  components  is 

90.91%  tungsten,  4%  iron,  3%  nickel,  and  2%  copper.  The  density  for  the  two 

3 3 

lots  was  17.0  g/cm  and  17.2  g/cm  . 

4.3  Preparation  Technique. 

The  first  step  in  preparing  a target  is  the  diffusion  welding  of  the 
junction.  The  pieces  a and  b (see  Fig.  4.2)  are  machined  to  the  following 
dimensions: 


(a)  63.5  mm  diam.  x 4 mm  thickness 

(b)  22.2  mm  diam.  x 6.4  mm  thickness. 

The  flat  surfaces  of  the  disks  are  then  lapped  flat  and  parallel  to  within  1 
micron  over  the  entire  area.  At  this  point,  the  surfaces  have  a damage  layer 
due  to  lapping  of  between  5 and  15  microns  thickness.  The  surfaces  to  be 
welded  are  now  carefully  polished  in  three  stages  on  a Buehler*  polishing 
wheel.  The  objective  in  the  first  stage  is  to  completely  remove  the  damage 
layer  due  to  lapping.  This  is  done  using  a nylon  polishing  cloth  loaded  with 
15  V diamond  polishing  paste  and  light  oil  designed  for  this  purpose.  The 
quality  of  the  results  depends,  to  a great  extent,  on  the  pressure  applied  to 
the  specimen,  and  on  the  amount  of  oil  and  diamond  paste  used.  At  the  com- 
pletion of  this  stage,  the  surface  should  appear  under  a 400x  microscope  as  a 
uniform  "scratched"  surface  where  only  the  damage  due  to  polishing  is  visible. 

In  the  second  polishing  stage  6 y diamond  polishing  paste  is  used  and 
the  finished  surface  should  again  have  a uniform  scratched  appearance,  but 
with  much  smaller  scratches.  Final  polishing  is  done  with  0.05  y alumina 


♦Buehler  Ltd.,  Evanston,  Illinois. 


polishing  powder  and  water  in  the  form  of  a slurry.  A napped  cloth  is  used 
in  this  stage  with  very  little  pressure.  After  this  final  polish,  the  sur- 
face under  a 400x  microscope  appears  free  of  scratches  and  polishing  damage. 
There  is  a slight  waviness  which  develops  on  the  surface  with  a wavelength 
of  10-50  y and  a depth  of  about  0.1  y to  1 y.  This  can  be  minimized  by 
reducing  the  pressure  and  time  of  the  final  stage  to  a minimum  necessary  to 
remove  the  scratched  layer.  This  technique  works  equally  well  for  copper  and 
constantan,  but  the  larger  surface  of  the  copper  makes  a great  deal  of  down- 
ward force  necessary  in  the  first  stage. 

The  next  step  in  the  welding  process  is  to  place  the  polished  sur- 
faces together  in  the  proper  position,  clamp  the  assembly  together,  and  heat 
it  in  a vacuum.  This  was  accomplished  here  by  using  two  flat  anvils  of  stain 
less  steel  with  eight  bolts  around  the  circumference.  Thin  sheets  of  mica 
are  placed  between  the  steel  anvils  and  the  specimen  to  keep  the  copper  or 
constantan  from  sticking  to  the  steel  anvils.  This  assembly  is  then  placed 
in  a vacuum  furnace,  and  heated  to  800°C.  It  takes  approximately  30  minutes 
to  reach  this  temperature,  and  then  the  temperature  was  maintained  for  60 
minutes.  The  furnace  was  then  turned  off  and  allowed  to  cool  for  about  four 

_5 

hours.  A vacuum  of  4 x 10  Torr  was  maintained  throughout  the  heating  cycle 

Quality  of  the  welded  interface  was  tested  in  several  ways.  The 
mechanical  strength  of  the  wel i in  a simple  tension  and  in  shear  appears  to 
be  about  the  same  as  the  bulk  copper  in  its  annealed  state.  No  elaborate 
mechanical  testing  was  done,  but  it  was  not  possible  with  brute  force  to 
separate  the  materials  even  with  severe  distortion  of  the  copper.  The  junc- 
tion was  carefully  sectioned  and  examined  under  a scanning  electron  micro- 
scope and  no  damage  or  impurities  were  detected. 

An  estimate  can  be  made  of  the  depth  of  interdiffusion  from  known 
diffusivities  of  nickel  in  copper.  Consider  the  situation  in  which  the 
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copper  and  constantan  are  separated  by  a plane  at  x = 0.  If  the  constantan 

39 

(45%  nickel)  is  to  the  left  of  the  dividing  plane.  Crank  gives  for  the 
concentration  of  nickel  for  x > 0 

C(x,t)  = i C erfc  — ^ . 

^ ° 2*^ 

Here  D is  the  diffusivity  of  nickel  in  copper  and  is  the  concentration 
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of  nickel  a,t  t = 0 and  x < 0.  dost  gives  the  value  for  D at  800°C  as 
-11  2 

4.84  X 10  cm  /sec.  Figure  4.3  is  a plot  of  C(x,  t = 1 hr)  for  a temper- 
ature of  800°C. 

After  the  welding  process  is  completed,  the  assembly  must  be  lapped 
square  again  due  to  a slight  deformation  of  the  copper.  The  piece  is  then 
carefully  aligned  in  a lathe,  and  the  constantan  is  machined  to  a dimension 
two  to  five  microns  larger  on  the  diameter  than  the  inside  dimension  of  the 
ceramic  tube.  This  dimension  is  nominally  15.25  mm,  which  allows  only  the 
central  portion  of  the  initial  weld  to  be  used,  thus  reducing  effects  due  to 
rounding  of  the  edges  of  the  constantan  during  polishing.  The  ceramic  and 
copper  guard  ring  is  assembled  by  a heat-shrinking  process.  The  inside  diam- 
eter of  the  copper  ring,  d,  is  machined  to  a dimension  two  to  five  microns 
smaller  than  the  outside  diameter  of  the  ceramic.  The  ceramic  tube  had  been 
previously  machined  round  to  a tolerance  of  less  than  2.5  microns.  The  copper 
is  then  heated  about  200°C  or  until  the  copper  and  ceramic  slip  together 
easily.  The  guard  ring  assembly  is  then  assembled  with  the  welded  copper- 
constantan  piece  by  cooling  the  copper-constantan  piece  in  liquid  nitrogen 
and  slipping  them  together.  Care  is  taken  throughout  this  process  to  main- 
tain the  parallelism  among  pieces.  The  entire  process  is  illustrated  in 
Fig.  4.2. 


4.4  Impact  Technique. 
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The  impact  experiments  were  done  using  the  four  inch  gas  gun  facility 
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at  Washington  State  University.  The  impactor  piece  was  mounted  on  the 
front  of  an  aluminum  projectile  as  illustrated  in  Fig.  4.4.  The  recess  which 
is  indicated  on  the  front  of  the  projectile  is  necessary  to  avoid  contact 
between  the  aluminum  projectile  face  and  the  copper  buffer  during  the  useful 
time  of  the  experiment.  For  projectile  velocities  greater  than  0.8  mm/usec  a 
nylon  and  syntactic  foam  projectile  was  used.*  The  design  of  this  projectile 
is  only  slightly  modified  from  that  reported  by  Peterson  et  ai.^^  Projectile 
velocity  is  measured  with  a set  of  four  electrically  charged  pins  which  are 
contacted  by  the  projectile  just  before  impact. 

The  target  assembly,  described  in  the  previous  section,  was  potted 
into  an  aluminum  target  ring  with  a calcium  carbonate  filled  casting  epoxy. 

A void  was  allowed  in  the  epoxy  around  the  soldered  oscilloscope  lead  on  the 
copper  buffer  to  avoid  a pressure  disturbance  in  this  area  during  the  experi- 
ment. A grounded  pin  and  an  electrically  charged  trigger  pin  contacted  the 
projectile  before  impact  and  were  used  to  start  the  oscilloscopes. 

Signals  in  the  one  to  twenty  millivolt  range  must  be  measured  with 
care  to  avoid  ground  loops  and  extraneous  noise.  The  signals  were  measured 
using  differential  amplifiers  in  a system  illustrated  in  Fig.  4.5.  This  type 
of  system  has  a 50  fl  input  impedance  in  each  leg  with  respect  to  ground,  thus 
giving  a 100  net  input  impedance.  The  analysis  of  the  thermoelectric  effect 
in  Chapter  3 assumes  an  open  circuit  or  infinite  impedance  measurement  of 
the  emf.  The  effect  of  this  difference  is  negligible  due  to  the  very  low 
internal  impedance  of  the  emf  source.  As  was  stated  in  Chapter  1,  the  emf 
source  here  can  be  treated  as  a voltage  source  with  an  associated  internal 

*Material  used  here  was  Scotchply  XP-241-34  Syntatic  Foam  supplied 
by  3M  Co.,  St.  Paul,  Minnesota. 


Velocity  Circuit 
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Tektronix 
7AI3  differential 
amplifier  in  a 
7904  oscilloscope 


Fig.  4.5. -“Schematic  Drawing  of  Recording  System. 


impedance.  The  ratio  of  the  measured  voltage  in  this  circuit  to  that  which 
would  be  measured  in  an  open  circuit  system  is 


1 


Input  Impedance 

Internal  Impedance  + Input  Impedance 


The  internal  impedance  of  this  thermoelectric  emf  source  was  measured  at  zero 
pressure  to  be  less  than  0.3  ohm.  This  would  lead  to  an  error  of  less  than 


0.3%. 


To  avoid  ground  loop  situations  in  which  the  shielding  of  a coaxial 
cable  is  carrying  a current,  the  shields  of  the  two  input  legs  are  connected 
together  and  to  the  scopes  only  and  are  isolated  from  all  other  possible 
ground  points.  It  is  not  unusual  in  this  facility  to  find  a difference  in 
potential  of  as  much  as  20  millivolts  between  points  that  are  nominally 
grounded.  Even  between  the  several  ground  points  on  a signal  oscilloscope  as 
much  as  a 10  millivolt  emf  can  be  measured.  To  contend  with  this  situation 


in  which  these  differences  in  ground  potential  produce  current  in  the  coaxial 
cable  shields,  ground  points  in  the  two  oscilloscopes  used  were  connected 
with  a heavy  copper  strap.  This  effectively  eliminated  the  problem. 

By  using  a differential  amplifier  to  measure  the  emf,  the  problem  of 
noise  on  the  scope  ground  is  minimized.  Not  only  is  there  a difference  in 
potential  between  various  ground  points,  but  there  is  also  as  much  as  five 
millivolts  of  60  cycle  noise  on  the  entire  chassis  ground  of  the  oscilloscope 
system.  Since  in  a differential  input  system  the  deflection  on  the  oscillo- 
scope is  proportional  to  the  difference  between  the  two  inputs  and  indepen- 
dent of  the  ground,  this  problem  can  be  completely  avoided.  Using  Tektronix 
model  7A13  amplifier  plug  in  units  in  model  7904  mainframes,  clean  signals 
could  be  measured  down  to  one  millivolt. 
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4.5  Comparison  of  the  Actual  Experiments  with  the  Idealized  Experiment. 

The  actual  configuration  described  in  section  4.1  and  the  idealized 
experiment  are  shown  for  comparison  in  Fig.  4.6.  The  similarities  of  these 
two  configurations  are: 

(1)  A plane  shock  wave  is  advancing  in  a cylinder  of  con- 
stantan  away  from  the  junction. 

(2)  The  copper-constantan  junction  is  in  a pressurized 
region  and  the  pressure  is  uniform  and  constant  over 
the  entire  junction. 

(3)  Pressure  relief  from  the  edges  of  the  constantan  behind 
the  shock  front  is  avoided.  (Assumed  condition  in 
idealized  experiment.) 

(4)  The  emf  is  measured  between  unshocked  regions  in  the 
copper  and  constantan. 

The  primary  differences  between  the  two  configurations  are: 

(1)  The  shock  wave  in  the  copper  is  not  a simple  plane  wave 
traveling  in  a cylinder,  but  is  rather  the  rounded  edge 
of  a plane  wave  traveling  radially  outward  in  the  copper 
disk. 

(2)  The  symmetry  with  respect  to  the  junction  is  not  main- 
tained in  the  real  experiment. 

(3)  The  pressure  in  the  copper  near  the  junction  is  attained 
by  a large  transmitted  shock  and  a small  reflected  shock 

tne  copper-constantan  interface.  The  reflected  shock 
»enr  wall  because  the  shock  impedances  of  copper  and 
— fwarir  equal  (see  Appendix  B). 

■•••  «•  •*9wmmn9'^  tr*  because  the  thermo- 
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Fig.  4.6. --Actual  Experimental  Configuration 
Compared  with  the  Idealized 
Configuration. 
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experiment  was  designed  to  preserve  the  ideal  conditions  in  the  constantan 
and  make  any  necessary  compromises  in  the  copper. 

The  experimental  problems  of  Chapter  3 were  avoided  in  these  experi- 
ments in  the  following  way: 

(1)  Electrical  noise  upon  circuit  closure  was  avoided  by 
having  the  recording  circuit  complete  before  the  shot. 

All  connections  are  either  soidered  or  welded  to  ensure 
good  contact. 

(2)  Demagnetization  effects  were  avoided  by  choosing  nonferro- 
magnetic materials.  Constantan  has  a relatively  high 
thermopower  of  about  50  yv/°C  at  ai d above  room  temperature 
and  is  nonferromagnetic.  The  primary  disadvantage  to  the 
use  of  this  material  arises  because  it  is  an  alloy.  Since 
thermopower  is  somewhat  affected  by  impurities  and  compo- 
sition, the  obvious  problem  of  reproducibility  among  samples 
of  different  manufacture  could  arise. 

(3)  A superheated  region  at  the  interface  was  avoided  or  mini- 
mized by  diffusion  welding  the  junction.  This  type  of 
junction  has  very  little,  if  any,  oorosity. 
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CHAPTER  5 

EXPERIMENTAL  RESULTS  AND  DISCUSSION 

Results  presented  here  are  divided  into  two  parts.  In  the  first  part 
the  guard  ring  was  modified  or  completely  eliminated.  This  allowed  for  pres- 
sure relief  radially  behind  the  shock  wave  in  the  constantan.  In  the  second 
part  the  results  of  a set  of  experiments  in  which  the  guard  ring  was  used  to 
maintain  one-dimensional  flow  behind  the  shock  will  be  presented.  Finally, 
our  conclusions  and  a proposal  for  further  study  will  be  given. 

5.1  Results  of  Experiments  without  Guard  Ring. 

Consider  again  the  idealized  experiment  of  Chapter  1.  If  the  edges 
of  the  rod  in  which  the  shock  wave  is  traveling  were  free  surfaces,  relief 
waves  would  advance  inward  from  the  edges  behind  the  shock  wave  to  reduce  the 
pressure.  It  was  assumed  in  the  analysis  that  no  pressure  relief  of  this 
type  was  allowed.  We  have  observed  that  this  restriction  has  a significant 
effect  on  the  measured  emf.  Four  experiments  were  conducted  in  which  the 
effects  of  this  restriction  can  be  seen.  The  experimental  parameters  and 
the  results  of  these  shots  (77-005,  77-026,  77-037,  and  77-074)  are  summarized 
in  Table  5.1,  and  Fig.  5.1.  A short  discussion  of  each  shot  will  follow. 

77-005  and  77-026.  These  experiments  were  conducted  as  described  in 
Chapter  4 except  that  no  guard  ring  was  used,  and  the  diameter  of  the  con- 
stantan was  19  mm  ± .1  mm.  The  second  shot  (77-026)  was  done  to  duplicate 
the  first,  and  was  identical  to  the  extent  possible.  Shock  pressure  was  cal- 
culated from  the  projectile  velocity  of  .7  mm/psec  to  be  144  i 4 kbars  in 
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both  shots.  This  would  give  a temperature  rise  in  the  shock  compression  of 
about  72°C  for  which  an  emf  of  about  3.5  mv  is  expected  as  calculated  in 
Chapter  1.  The  voltage-time  profiles  (Fig.  5.1a  and  b)  are  characterized  by 
a negative  step  at  the  instant  the  shock  wave  reaches  the  copper-constantan 
interface  of  -3  ± .5  mv  and  -4  ± .2  mv,  respectively.  This  step  is  followed 
by  a negative  ramp.  In  both  cases  the  slope  of  the  ramp  decreases  with  time, 
but  initial  slopes  are  different  (-24  mv/ysec  and  -14  mv/ysec,  respectively). 
The  ramp  is  followed  by  a positive  step  of  2.5  mv  and  3 mv,  respectively, 
when  the  shock  wave  reaches  the  back  surface  of  the  constantan. 


TABLE  5.1 .--Results  of  Experiments  with  Radial  Pressure  Relief 


Shot 

Number 

Junction 

Materials 

Projectile 

Velocity 

(mm/ysec) 

Shock 

Pressure 

(kbar) 

Initial  Step 
Ampl itude 
(mv) 

Initial 
Slope 
of  Ramp 
(mv/ysec) 

77-005 

OFHC  Copper, 
A-Constantan 

0.70  ± 

.01 

144  ± 4 

-3  ± .5 

-24 

77-026 

OFHC  Copper, 
A-Constantan 

0.70  ± 

.01 

144  ± 4 

-4  ± .2 

-14 

77-037 

OFHC  Copper, 
A-Constantan 

1.31  ± 

.02 

300  ± 8 

-15+2* 

-75 

77-074 

OFHC  Copper, 
B-Constantan 

0.70  ± 

.01 

144  ± 4 

-3.2  ± .2 

-13 

*The  step  height 

is  difficult 

to  determine 

due  to  the 

steepness  of 

the  following  ramp. 

77-037.  The  projectile  velocity  was  increased  in  this  experiment  and 
was  measured  to  be  1.31  ± .02  mm/ysec.  The  dimensions  of  the  target  and 
impactor  were  scaled  down  to  facilitate  the  high  projectile  velocity.  The 
constantan  diameter  was  the  same  at  19  mm,  the  impactor  diameter  was  25  mm, 
the  copper  buffer  diameter  was  35  mm,  and  the  copper  buffer  thickness  was 


1.7  mm  (nominal  dimensions  with  .1  mm  tolerance).  The  pressure  was  calcu- 
lated to  be  300  ± 8 kbars  which  would  give  a temperature  rise  in  the  shock 
compression  of  about  194°C  and  an  estimated  thermal  emf  of  about  9.5  mv.  The 
voltage-time  profile  (Fig.  5.1c)  is  quite  different  for  this  shot  than  for 
the  previous  two  shots.  A step  in  emf  is  not  clearly  distinguishable  from 
the  following  ramp.  In  the  original  data  trace  a slight  change  in  slope  can 
be  seen  at  about  -15  mv.  The  magnitude  of  the  emf  does  not  increase  monoton- 
ically  during  the  shock  transit  through  the  constantan  as  in  shots  77-005  and 
77-026,  but  reaches  a maximum  of  -36.5  mv  about  .2  ysec  before  the  shock  wave 
reaches  the  back  surface  of  the  constantan. 

77-074.  In  this  shot  the  guard  ring  of  ceramic  and  copper  described 
in  Chapter  4 was  used,  but  unlike  the  configuration  in  Fig.  4.1  the  con- 
stantan thickness  (B)  was  13  ± .03  mm  and  the  guard  ring  thickness  (E)  was 
only  5.16  ± .03  mm.  We  had  tentatively  determined  at  this  time  that  the 
ramping  in  shots  77-005,  77-026,  and  77-037  was  related  to  the  violation  of 
the  one-dimensional  strain  restriction  of  the  idealized  experiment.  This  was 
supported  by  the  steady  emf  (no  ramping)  observed  in  another  shot  77-064  in 
which  the  guard  ring  was  used  as  described  in  Chapter  4 (see  section  5.2). 

The  specific  purpose  of  77-074  was  to  determine  whether  the  pressure  relief 
at  the  junction  itself  causes  the  ramping,  or  if  it  could  be  observed  when 
the  radial  flow  was  allowed  just  behind  the  shock  front,  but  not  at  the  junc- 
tion. The  projectile  velocity  and  shock  pressure  were  the  same  as  for  shots 
77-005,  and  77-026.  The  resulting  voltage-time  profile  (Fig.  5. Id)  shows  an 
initial  step  of  -3.2  ± .2  mv  at  the  instant  the  shock  wave  reaches  the  copper- 
constantan  interface  and  a steady  emf  until  the  shock  wave  reaches  the  back 
surface  of  the  guard  ring.  At  that  point  a ramping  of  -13  mv/usec  is  obser- 
ved. This  indicates  that  the  ramping  is  associated  with  the  relief  process 
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at  the  edges  in  the  constantan,  and  is  not  due  to  the  interaction  of  the 
relief  waves  with  the  copper-constantan  interface. 

These  results  indicate  that  in  the  case  where  simple  uniaxial  strain 
conditions  exist,  the  net  emf  around  the  thermoelectric  circuit  is  constant 
as  the  shock  wave  progresses  through  the  constantan.  In  the  case  where  two- 
dimensional  flow  exists  behind  the  shock  front  in  constantan  due  to  relief  at 
the  edges,  the  emf  around  the  circuit  is  greater  in  magnitude  and  unsteady, 
at  least  under  these  experimental  conditions.  In  the  experiments  where  no 
guard  ring  was  present,  the  emf  measured  when  the  shock  first  pressurizes  the 
junction,  and  before  relief  waves  from  the  edges  are  established,  compares 
directly  with  the  steady  emf  measured  when  a guard  ring  was  used,  as  shown  in 
Table  5.1. 

We  propose  three  possible  causes  for  this  difference  in  electrical 
response  between  these  two  types  of  experiments. 

(1)  High  velocity  radial  flow  due  to  radial  pressure  relief 
which  would  tend  to  rapidly  change  the  self-inductance 
of  the  current-carrying  conductor.  The  entire  circuit 
then,  when  taken  as  a unit,  would  have  an  unsteady  cur- 
rent flow  even  for  a constant  emf  source. 

(2)  Relativistic  effects  usually  neglected  in  a moving  con- 
ductor. This  is  unlikely  because  of  the  small  speeds 
involved,  but  is  included  for  completeness. 

(3)  Some  unknown  transport  property  of  the  electron- 
lattice  system  making  up  the  conductor  when  two- 
dimensional  flow  is  present  behind  the  shock  front. 

Calculations  have  been  done  In  Appendix  F which  Indicate  that  the 
effects  Mentioned  in  the  first  two  possible  explanations  are  far  too  .xwli 
ween  ( a<  ■ •wn*  -Wtservee 


I 
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A serious  study  of  the  possibilities  involved  in  the  third  suggested 
explanation  has  not  yet  been  undertaken.  It  should  be  noted,  however,  that 
the  absolute  thermopower  is  affected  to  varying  degrees  by  both  hydrostatic 
pressure  and  one-dimensional  strain  in  static  situations.  Both  of  these 
effects  can  be  understood  at  least  qualitatively  on  the  basis  of  modern  trans- 
port theory  as  it  applies  to  thermoelectric  phenomena.  We  feel  that  further 
study  of  this  problem  is  warranted  on  the  basis  of  these  results. 

5.2  Results  of  Experiments  with  Guard  Ring. 

Six  additional  measurements  were  made  with  the  experimental  design  of 
Fig.  4.1  and  4.4,  including  the  ceramic  and  copper  guard  ring.  The  only 
parameters  that  were  varied  in  these  experiments  were  the  shock  pressure  and 
the  constantan  material.*  The  results  of  one  additional  shot,  77-074,  men- 
tioned in  section  5.1,  in  which  the  constantan  was  longer  than  the  guard  ring, 
compare  directly  for  times  less  than  the  shock  transit  time  through  the  guard 
ring,  and  will  be  included  in  these  results. 

Results  of  these  experiments  are  listed  in  Table  5.2,  The  pressure 
behind  the  shock  was  calculated  for  each  of  the  shots  using  the  Rankine- 
Hugoniot  jump  conditions  and  empirical  Ug-Up  relations  for  both  copper  and 
constantan  (see  Appendix  B).  The  uncertainties  indicated  for  the  measure- 
ments are  based  on  the  quality  of  the  record  and  its  character;  i.e.,  where 
records  are  noisy,  or  the  emf  is  not  steady  but  varies  slightly  with  time, 
the  uncertainty  attributed  is  larger.  The  uncertainties  in  each  case  are 
larger  than  could  be  attributed  to  the  recording  system  or  data  handling 
system;  therefore,  the  uncertainty  listed  reflects  only  the  character  of  the 
observed  signals  rather  than  our  ability  to  record  the  signal. 


*In  shot  77-OM  a disk  of  copter  aas 
tht  re*staMta«<  as  a«*  eioctrode  twiteod  sv 
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TABLE  5.2. --Results  of  Experiments  with  Uniaxial  Strain 
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The  voltage-time  profiles  for  all  of  these  experiments  are  given  in 
Figs.  5.2  and  5.3.  The  initial  negative  step  in  emf  in  each  of  the  records 
takes  place  when  the  shock  wave  reaches  the  copper-constantan  junction.  In 
the  shots  done  near  150  kbars,  the  emf  remains  approximately  constant  until 
the  shock  wave  reaches  the  back  free  surface  of  the  constantan.  The  rise 
times  for  the  initial  voltage  step  are  included  in  Table  5.2  for  each  shot. 

The  rise  times  vary  considerably,  but  the  short  rise  time  of  7 ns  observed  in 
shot  77-076  indicates  that  the  response  time  of  the  thermoelectric  circuit  is 
short  compared  with  the  typical  value  observed  for  rise  time.  A great  deal 
of  care  is  taken  in  an  effort  to  maintain  parallelism  between  the  advancing 
plane  shock  and  the  junction  plane.  However,  neither  the  tilt  of  the  projec- 
tile with  respect  to  the  target  nor  the  actual  parallelism  between  target  face 
and  junction  were  measurable  in  these  experiments. 

In  the  two  highest  pressure  shots,  77-079  and  77-092,  the  emf  became 
unsteady  and  increasing  in  magnitude  after  about  300  ns.  This  has  been 
attributed  to  the  mechanical  behavior  of  the  ceramic  insulator  which  makes  it 
a poorer  impedance  match  to  the  copper  and  constantan  at  high  pressures  than 
at  low  pressures  as  discussed  in  section  4.2.  The  lower  pressure  in  the 
ceramic  allows  pressure  relief  at  the  edges  of  the  constantan  and  produces  an 
effect  similar  to  that  seen  in  the  unsupported  shots. 

5.3  Comparison  of  Results  with  Theory  and  Discussion. 

The  emfs  observed  for  the  experiments  discussed  in  section  5.2  can 
now  be  compared  directly  with  calculations  based  on  the  thermoelectric  theory, 
known  aaiplrlcal  relations  for  thermopower,  and  temperature  estimates  discus 
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2.— Voltage-Time  Profiles  for  Shots  77-064,  77-076,  and  77-079 
The  large  tick  marks  indicate  the  time  at  which  the  shock 
wave  reached  the  hack  surface  of  the  constantan. 
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Voltaqe-TIme  Profiles  for  Shots  77-084,  77-085,  and  77-092 
The  1arqc  tick  marks  Indicate  the  tine  at  which  the  shock 
wave  reached  the  back  surface  of  the  constantan 


one-dimensional  loading  of  the  junction.  Therefore,  the  initial  steps  for 
those  experiments  can  be  compared  to  the  estimates  of  Chapter  1 as  well. 
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In  Fig.  5.4,  the  data  for  shots  in  which  A-constantan  was  used  are 
shown  graphically.  For  shots  in  which  the  guard  ring  was  omitted  the  data 
point  represents  only  the  initial  step.  The  lower  curve  represents  an  evalu- 
ation of  the  formula  in  Eq.  (1.43)  for  the  temperature  estimates  of  section 
1.3.  The  upper  curve  represents  a similar  evaluation  where  the  pressure  cor- 
rection to  the  thermopower  is  neglected. 

The  B-constantan  was  used  primarily  in  the  development  stage  of  the 
experimentation,  and  no  attempt  was  made  to  vary  the  pressure  when  it  was 
used.  The  results  for  the  three  shots  reported  in  sections  5.1  and  5.2  for 
this  type  of  material  can  be  compared  with  a calculated  emf  value  based  on 
the  considerations  of  Chapter  1 and  the  standard  thermopower  of  constantan.*^^ 
The  initial  emf  step  as  well  as  the  calculated  emf  are  given  in  Table  5.3 
for  the  shots  involving  B-constantan. 

The  basic  correlation  between  the  measured  emf  and  the  emf  calculated 
from  shock  compression  temperature  rise  indicates  that  the  observed  emf  is 
indeed  of  thermoelectric  origin.  In  this  context,  this  co. relation  can  be 
evaluated  more  easily  by  comparing  the  calculated  temperature  rise  in  the 
shocked  region  with  the  "measured  temperature"  which  will  be  defined  as  the 
temperature  rise  necessary  to  ^tve  the  observed  results  if  the  junction  acts 
a^  a thermocouple.  In  Table  5.3,  the  "measured  temperature."  based  on  a pres- 
sure "^nrrected  theimuuoMer  a<.  dHcussed  ’r  sec*>oe  ’ is  :i»er  *nr  *■ 

^ V ' * ^ • r # s ' l ***•  >1^ 


77-037, 


TABLE  5.3. --Comparison  of  Calculated  and  Measured  Temperatures 
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temperature  and  the  "measured  temperature"  is  22%.*  Three  possible  sources 
of  error  involved  in  this  comparison  are  the  error  in  the  measured  emf,  the 
uncertainty  in  the  calculated  shock  compression  temperature,  and  the  validity 
of  the  pressure  correction  to  the  thermopower.  Only  the  errors  involved  in 
the  emf  measurement  for  each  shot  are  indicated  explicitly  as  an  uncertainty 
in  the  "measured  temperature."  The  shock  compression  temperature  calculation 
using  the  assumptions  that  y/v  and  are  constant  is  believed  to  be  cor- 
rect within  about  10%.  The  temperature  in  the  shocked  region  could  be  as 
much  as  5%  higher  than  the  value  given  due  to  plastic  work  if  the  yield 
stress  were  significantly  larger  than  the  value  obtained  from  quasistatic 
measurements.  The  pressure  correction  to  the  thermopower  is  extrapolated 
from  hydrostatic  thermopower  corrections  measured  to  72  kbars.  The  actual 
magnitude  of  this  correction  at  any  given  pressure  would  depend  somewhat  on 
the  pressure  dependence  of  temperature  in  the  shock  front  (see  section  1.2). 
However,  since  this  correction  is  based  theoretically  on  the  effects  of  the 
change  in  atomic  volume,  the  fact  that  the  pressure  is  reached  through  a 
shock  mechanism  rather  than  hydrostatically  should  not  make  the  correction 
completely  invalid.  At  150  kbars  this  correction  increases  the  "measured 
temperature"  by  only  5%. 

Considering  the  10%  uncertainty  in  the  calculated  temperature,  the 
measured  emf  in  the  experiments  corresponds  to  a temperature  that  may  be 
within  10%  of  the  actual  shock  compression  temperature.  At  least  three  pos- 
sible explanations  exist  with  respect  to  the  remaining  discrepancy.  First, 
it  is  possible  that  the  calculated  temperature  is  10-20%  lower  than  the 
actual  value  and  that  the  "measured  temperature"  is  more  accurate  than  the 
calculated  value.  Attempts  have  been  made  to  determine  the  shock  compression 

*This  is  the  RMS  value  neglecting  shot  77-037. 
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temperature  of  opaque  materials  by  measuring  the  residual  temperature  after 

44-46 

expansion  to  zero  pressure.  The  techniques  used  are  based  on  measuring 

the  surface  radiation  temperature  or  resistivity  of  the  unloaded  material. 
Sufficient  data  are  not  available  at  this  time  to  adequately  examine  this 
possibility.  Secondly,  the  10-20%  discrepancy  may  be  due  to  a slight  change 
in  the  thermopower  in  shock  compression.  This  would  represent  a change  in 
the  heat  transport  properties  of  the  conduction  electrons,  and/or  the  crystal 
lattice  from  the  static  situation.  Finally,  the  temperature  in  the  immediate 
vicinity  of  the  junction  interface  may  be  somewhat  higher  than  the  bulk  tem- 
perature. We  have  shown  in  section  3.3  that  the  temperature  measured  in  this 
type  of  thermoelectric  circuit  may  be  that  which  exists  within  a few  microns 
of  the  interface  rather  than  in  the  bulk  material  behind  the  shock. 

It  is  our  belief  that  the  interface  temperature  could  easily  be  10-20% 
higher  than  the  bulk  temperature,  and  that  this  is  the  most  probable  expla- 
nation of  the  discrepancy  between  the  calculated  and  measured  values  of 
temperature.  The  diffusion  welding  process  used  here  reduces  the  difference 
between  a real  junction  and  an  ideal  junction,  but  some  imperfections  undoubt- 
edly exist  which  could  lead  to  a higher  temperature  in  this  region.  In  this 
context,  the  diffusion  welded  junction  is  one  step  closer  to  the  ideal  junc- 
tion in  which  a plane  divides  the  two  materials. 

5.4  Conclusions  and  Recommendations. 

5.4.1  Conclusions 

The  initial  objective  to  determine  whether  the  previously  reported 
anomalous  and  generally  nonreproducible  results  were  due  to  experimental 
artifacts  or  fundamental  principles  has  been  achieved.  Our  results  indicate 
that  any  fundamental  anomaly  which  exists  for  one  dimensional  strain  is  less 
than  10-20%  of  the  observed  signals  in  the  piressure  range  from  140-360  kbars. 
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In  a study  of  the  effects  of  an  imperfect  interface  layer  between  the  metals 
due  to  surface  preparation,  we  have  shown  that  the  measured  emf  in  a thermo- 
electric circuit  is  sensitive  to  the  junction  interface  temperature,  which 
might  be  several  times  greater  than  the  shock  compression  temperature  of  the 
bulk  material.  We  therefore  conclude  that  the  large  emf  values  previously 
I I reported,  as  well  as  the  general  nonreproducibility  among  experiments,  were 

I due  primarily  to  the  effects  of  poor  junction  interface  conditions.  Using  a 

’ model  of  a porous  layer  between  the  materials,  it  was  shown  that  the  temper- 

I : ature  in  this  region  is  nearly  proportional  to  the  pressure  in  a single  shock 

I compression.  This  may  explain  the  linear  emf  pressure  relations  previously 

' 19  28 

observed.  In  studies  by  Buzhinskii  et  ^. , and  Nesterenko,  in  which 

particular  attention  was  given  to  surface  preparation,  some  of  the  resulting 

emf  values  are  within  20%  of  the  predicted  values  and  although  a great  deal 

I of  scatter  in  the  data  is  present  in  both  studies,  this  evidence  supports  the 

t 

I conclusion  that  no  fundamental  anomaly  exists  and  that  high  emf  values  are 

( 

I due  to  poor  surface  conditions. 

I We  have  shown  that  in  experiments  designed  to  measure  the  response  of 

I a bimetallic  junction,  it  is  possible  to  inadvertently  measure  an  induced  emf 

^ caused  by  shock  demagnetization  when  a ferromagnetic  material  such  as  nickel 

I is  used  as  one  of  the  elements.  This  may  be  responsible  for  some  of  the 

I nonreproducibility  evident  in  previous  investigations  in  which  ferromagnetic 

I materials  were  used. 

I ' In  experiments  in  which  radial  pressure  relief  was  allowed  behind  the 

I J shock  front,  the  emf  was  observed  to  be  unsteady  and  increasing  with  time. 

t I 

I j No  adequate  explanation  of  this  observation  has  yet  been  found.  We  have 

I ■ 

I shown,  however,  that  the  unsteady  emf  is  directly  related  to  the  two- 

f I dimensional  flow  caused  by  radial  pressure  relief  behind  the  shock.  This 
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effect,  having  now  been  isolated,  may  be  considered  as  a possible  source  of 

the  nonreproducibility  present  in  some  previous  experiments.  In  particular, 

20 

the  relatively  small  emf  followed  by  a ramping  observed  by  Lascar  and  Duage 
could  be  explained  as  a "normal"  thermoelectric  response  coupled  with  this 
two-dimensional  effect. 

[ ig 

I The  one  oscillogram  reported  by  Buzhinskii  et  does  not  indicate 

i a ramping  effect  similar  to  that  observed  here.  The  same  basic  experimental 

geometry  was  used  in  that  study  as  was  used  for  the  experiments  reported 
here;  however,  in  some  experiments  the  nickel  was  in  the  shape  of  a 120“  cone 
which  would  modify  the  relief  wave  situation.  The  authors  do  not  indicate 
which  experiments  were  done  in  this  way  or  for  which  experiment  the  oscillo- 
graph was  observed.  Since  the  surfaces  and  geometry  were  varied  and  the 
data  have  a great  deal  of  scatter,  it  is  impossible  at  this  point  to  determine 
how  these  results  compare  with  those  observed  here. 

Finally,  we  have  shown  that  the  voltage  response  of  a carefully  pre- 
pared and  diffusion-welded  copper-constantan  junction  to  one-dimensional 
shock  compression  over  a range  in  pressure  from  145  kbars  to  360  kbars  is 
within  10-20%  of  the  predicted  value.  This  prediction  is  based  on  a pressure 
corrected  thermoelectric  response  to  the  shock  compression  temperature  of  the 
bulk  materials.  Further  study  is  necessary  to  determine  conclusively  whether 
this  remaining  discrepancy  is  due  to  the  inaccuracy  of  the  temperature  calcu- 
lation or  of  the  temperature  measurement.  At  this  point  we  believe  that  the 

t 

S discrepancy  is  more  probably  due  to  a temperature  at  the  junction  interface 

it 

I 

I which  is  higher  than  that  of  the  bulk  materials  due  to  imperfections  in  the 

r 

welded  interface. 

t. 

i 

\ 

i 

\ 

I 

i 

[ 

1 

^ — - - - 


i 
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• 5.4.2  Recommendations 

I 

i We  recommend  that  further  study  be  done  to  determine  the  cause  of  the 

I 

; decreasing  unsteady  emf  observed  when  radial  pressure  relief  is  allowed 

behind  the  shock  front.  If  this  effect  is  indeed  due  to  a fundamental  change 
in  the  transport  properties  of  the  material  under  these  conditions,  charac- 
I terization  of  this  phenomenon  could  considerably  advance  basic  understanding 

in  this  area. 

Further  study  is  also  recommended  to  determine  the  cause  of  the 
remaining  discrepancy  between  predicted  and  measured  temperatures.  Since  the 
measured  temperature  is  that  which  exists  in  the  immediate  vicinity  of  the 
junction  interface,  the  possibility  of  a higher  temperature  due  to  defects  at 

1 

1 

I or  near  the  weld  should  be  eliminated.  The  use  of  a smaller  junction  area  is 

recommended  to  facilitate  the  diffusion  welding  process. 

The  use  of  a highly  compressible  metal  as  one  element  in  this  junction 
would  generate  a much  higher  expected  temperature  and  make  the  imperfections 
i at  the  weld  negligible.  The  discrepancy  between  the  measured  and  predicted 

temperatures  should  then  be  smaller  if  this  discrepancy  is  due  only  to  defects 
I at  the  interface. 

I 

r Finally,  we  recommend  that  the  feasibility  be  studied  of  using  a 

i 

I thermocouple  as  a temperature  transducer  in  shock  compressed  materials  other 

! than  the  thermocouple  elements.  We  have  shown  that  the  thermoelectric  effect 

f 

i is  not  altered  to  any  large  degree  by  a one-dimensional  shock  compression  and 

i therefore  might  be  employed  in  a shock  environment  to  determine  the  tempera- 

ture of  the  shocked  state  in  surrounding  material. 


j 

f 

i 

k 
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APPENDIX  A 

MEASUREMENT  OF  THE  THERMOELECTRIC  POWER  OF  AN  OFHC  COPPER 
AND  A-CONSTANTAN  THERMOCOUPLE 

The  thermoelectric  power  for  the  materials  OFHC  copper 

and  A-constantan  was  determined  by  measuring  the  potential  difference 
as  in  Fig.  1.5.  The  data  taken  for  various  values  of  T are  given  in 
Table  A.l  where  T^  = 23  ± 2°C.  The  standard  thermocouples  were  made  from 
wire  manufactured  by  Leeds-Northrop  Co.  with  a rated  accuracy  of  ± 3/4%. 

TABLE  A.l .--Thermopower  Data  at  Zero  Pressure 


(T-Tq) 

(°C) 

<J>E-F 

(mv) 

77 

3.80 

85 

4.24 

111 

5.63 

149 

7.77 

154 

8.20 

174 

9.27 

191 

10.30 

199 

10.76 

248 

14.00 

276 

15.85 

289 

16.54 

310 

17.84 

Using  the  method  of  least  squares  these  data  were  fit  to  the  expression 
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♦e-f(T-T„)  = Q(T-T„)  + R(T-T„)2 

where 

Q = 0.0475  mv/°C 
R = 3.4  X 10“^  mv/°C^  . 

The  standard  error  of  estimate  is  0.095  mv.  In  Fig.  A.l  the  data  taken  for 
<I>£_p(T-Tq)  and  the  associated  least  squares  fit  are  given. 


The  solid  line  represents  a least  squares  fit  to  the  data. 
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The  temperature  behind  a single  shock  can  be  calculated  as  outlined 
in  section  3.3.  The  temperature  as  a function  of  final  shock  pressure  was 
found  by  numerically  integrating  the  expression 


j T = e-'^(''-''o)'^o  + e-^''  j e'^'''F(v' )dv- 

(where 

^ = ^o/'o  ’ 


I 


F(v)  = 


■ y dv^  2C  " ^x  " 


- V dp 

r — (-  -sr>i 

X HEL 


T = {a  - a )/2  (maximum  resolved  shear  stress) 

A y 

Here  y is  the  yield  stress  in  uniaxial  stress,  y is  the  shear  modulus, 
is  the  component  of  pressure  in  the  shock  propagation  direction  and  is 
the  specific  heat  at  constant  volume.  The  subscript  HEL  refers  to  the  value 
of  the  parameter  at  the  Hugoniot  elastic  limit. 

Work  hardening  can  be  included  by  letting  t be  a function  of 
volume.  In  simple  uniaxial  compressive  stress  the  yield  stress  can  be 
expressed  as 
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Y = Yp  + - e°)  (B.2) 

where  is  the  natural  strain  in  uniaxial  stress  and  H is  a constant. 

47 

Fowles  has  shown  that 

where  is  the  natural  strain  in  uniaxial  strain  and  K is  the  isothermal 
bulk  modulus.  If  the  variation  in  K is  neglected  we  have 

S " f (s  - 

Equation  B.2  can  now  be  expressed  as 

(Y  - Yq)/2  = [|  (e^  - e°)]/(l  - . (B.5) 

U 

For  reasonable  values  of  H and  K the  term  ^ is  negligible  compared  to 
1.  The  natural  strain  for  the  one-dimensional  strain  case  can  be 
expressed  as  a function  of  volume  as 

^ ''o 

Using  the  von  Mises  yield  criteria  for  uniaxial  strain, 


Eg.  (B.5)  becomes 


3 


! 

i 


(B.7)  I 


To  calculate  the  temperature  using  Eq.  (B.l)  the  pressure  must  be 
known  as  a function  of  volume  along  the  Hugoniot.  As  the  integration  pro- 
ceeds the  values  of  p and  dp/dv  are  calculated  using  a - Up  relation 
and  the  Rankine-Hugoniot  relations.  Using  the  relation 


The  Ug  - Up  relations  used  for  these  calculations  and  to  determine  the 
final  pressure  in  a given  experiment  are: 


OFHC  Copper^® 

Cq  = 3.94  mm/ysec 

s 

= 1.489 

Tungsten  Alloy 

Cq  = 3.97  mm/ysec 

s 

= 1.11 

A-Constantan 

Cq  = 4.24  mm/ysec 

s 

= 1.55 

B-Constantan 

Cq  = 4.20  mm/ysec 

s 

- 1 .55 

The  relation  for  the  tungsten  alloy  was  obtained  by  measuring  the  shock  speed 
in  two  symmetric  impact  experiments  at  about  240  kbars  and  355  kbars.  Also 
the  zero  pressure  intercept  was  estimated  from  ultrasonic  sound  speed  measure- 
ments. The  relation  for  A-constantan  was  obtained  by  measuring  the  shock 
speed  in  shots  77-005,  77-026,  77-037,  77-079,  and  77-084.  From  the  data 
trace  (see  Chapter  5)  one  can  obtain  the  shock  transit  time  through  the  con- 
stantan.  With  this  and  the  known  Hugoniot  of  copper  one  can  determine  the 
Us  - Up  point  for  each  experiment.  The  standard  error  of  estimate  for  these 
five  shots  using  the  least  squares  fit  given  above  is  0.05  mm/ysec.  The 
relation  for  B-constantan  was  obtained  by  shifting  the  curve  for  A-constantan 


to  intercept  the  zero  pressure  bulk  sound  speed  obtained  from  ultrasonic 
sound  speed  measurements. 
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The  values  for  the  remaining  parameters  necessary  in  this  calculation 
which  were  used  here  are  given  in  Table  B.l.  Temperatures  calculated  using 
Eq.  (B.l)  are  given  in  Table  B.2.  For  comparison,  temperatures  are  listed 
with  (A  columns)  and  without  (B  columns)  the  work  hardening  correction.  In 
the  last  column  the  calculated  values  of  Rice  e^al^.^*^  for  copper  are 
included. 


TABLE  B.l .--Constants  Used  in  Temperature  Calculations 


Parameter 

OFHC  Copper 

A-constantan 

B-constantan 

★ 

Y 

1.99 

2.26 

2.06 

''hEL^'o 

0.99974 

0.999382 

0.99187 

Yq  (kbar)^ 

-.2 

-1 

-1 

Tq  (kbars)' 

0.125 

0.5 

0.5** 

H (kbars)^ 

15.5 

3.74 

3.74** 

C (kbars)* 

480. 

810. 

615. 

pQ  (g/cm^) 

8.93 

8.93 

8.82 

(10^  ergs/gr°K)* 

3.73 

3.82 

3.83 

V* 

.35 

.30 

.32 

ultrasonic  longitudinal 
sound  speed  (mm/ysec) 

4.81 

5.60 

5.187 

ultrasonic  transverse 
sound  speed  (mm/ysec) 

2.32 

3.01 

2.64 

*Derived  from  ultrasonic  sound  speed  measurements  done  here  and 

handbook  values  of  C„  in  the  case  of  C . 

P V 

^From  quasi  static  yield  stress  measurements  done  here. 

**Value  for  A-constantan  was  used. 


^ 


*The  "A"  columns  are  values  with  work  hardening, 

i **The  "B"  columns  are  values  without  work  hardening. 

I + 

» • 1 r\ 
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APPENDIX  C 

ELECTRICAL  NOISE  AT  IMPACT 

Four  experiments  were  done  specifically  to  investigate  the  problem  of 
electrical  signal,  produced  when  a thermoelectric  circuit  is  completed  by  the 
‘ impact.  The  experimental  configuration  is  illustrated  in  Fig,  C.l,  and  the 

voltage  time  profile  for  each  shot  is  given  in  Fig.  C.2.  Since  each  experi- 
ment was  somewhat  different  a short  description  of  each  will  follow, 

76-073.  The  configuration  is  symmetric  with  respect  to  the  impact 
plane  except  for  the  Nichrome  V center  element  in  the  target.  This  material 
' was  used  here  to  simulate  the  constantan  since  it  has  similar  mechanical  and 

I magnetic  properties  and  is  more  readily  available.  The  impactor  assembly  was 

■ mounted  on  an  aluminum  projectile  (see  section  4.4).  The  oscilloscopes  were 

triggered  about  10  usee  before  impact.  The  emf  was  measured  between  the 
inner  and  outer  pieces  of  the  target  assembly  with  a differential  measuring 
system  as  described  in  section  4.4.  The  magnetic  field  in  the  impact  area 
was  about  0.15  Gauss  in  a direction  perpendicular  to  the  direction  of  motion. 

The  projectile  velocity  for  this  series  of  shots  was  0.6  mm/usec  ± 0.02  mm/ysec. 

76-075.  In  this  experiment  the  magnetic  field  in  the  area  of  impact 
was  lowered  to  less  than  0.005  Gauss  in  any  direction  with  Helmholtz  coils. 

In  all  other  respects  the  experiment  was  the  same  as  76-073. 

76-078.  The  projectile  in  this  shot  was  modified  by  placing  a one 


inch  plastic  piece  on  the  front  of  the  projectile  between  the  aluminum  body 
and  the  copper  active  elements.  This  was  done  to  reduce  the  effects  of  a 
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elements 


Fig.  C.l .--Concentric  Ring  Experimental  Configuration. 
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moving  conductor  in  a magnetic  field  near  the  impact.  In  all  other  respects 
the  experiment  was  the  same  as  76-073. 
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76-086.  The  center  element  of  this  target  was  made  of  copper  which 
makes  the  projectile  and  target  assemblies  indentical.  Again,  the  experiment 
was  the  same  in  all  other  respects  as  76-073. 

The  objective  of  this  study  was  to  determine  the  cause  of  this  type 
of  noise  and  eliminate  it  if  possible.  By  lowering  the  magnetic  field  in  the 
impact  area  the  amplitude  of  the  signal  was  reduced  as  seen  in  76-075.  We 
concluded  that  the  pre-impact  signal  and  associated  noise  at  and  just  after 
impact  were  probably  due  to  an  interaction  of  the  moving  conductors  with  the 
magnetic  field.  We  determined  that  if  the  remaining  signal  observed  in  76-078, 
where  the  field  was  reduced  to  0.005  Gauss,  is  due  to  the  residual  magnetic 
field,  it  would  be  impossible  to  eliminate  this  type  of  noise  in  our  facility. 
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; APPENDIX  D 

1 

\ SHOCK  DEMAGNETIZATION 

I 

I Eight  experiments  were  conducted  in  a study  of  the  effects  of  shock 

^ demagnetization  in  bimetallic  junction  experiments.  The  basic  experimental 

I technique,  based  on  the  use  of  a four-inch  light  gas  gun,  is  detailed  in 

S section  4.4.  Instrumentation  involved  a combination  of  oscilloscopes  used  in 

: 

i 

■ a differential  mode  as  indicated  in  the  schematic  drawing  of  Fig.  4.5.  The 

1 

I material  used,  except  where  noted,  was  nickel-200,  formerly  called  "A"  nickel, 

I which  is  nominally  99.4%  nickel.  In  all  experiments,  except  76-054  and 

I 76-051,  the  magnetization  state  was  not  altered  from  the  supplied  condition 

j and  was  assumed  to  be  random  in  direction  with  respect  to  the  axis  of  shock 

I propagation.  Since  each  of  the  experiments  was  performed  with  a different 

objective  in  view,  a short  description  of  each  will  follow,  with  the  results 
summarized  in  Table  D.l. 

76-032.  A concentric  guard  ring  geometry  was  used  as  illustrated  in 
‘ Fig.  C.l.  The  insulator  used  was  machined  from  solid  teflon.  The  outer  ring 

used  in  this  case  for  support  and  to  complete  the  circuit,  was  made  from 
OFHC  copper.  The  configuration  is  symmetric  with  nickel  impacting  nickel  as 
active  elements  so  that  no  thermoelectric  emf  is  expected.  However,  a non- 
zero signal  was  observed.  The  voltage  versus  time  profile,  as  seen  in 
Fig.  D.l,  indicates  a zero  voltage  at  impact  with  a continuously  increasing 
magnitude  of  negative  sign.  This  result  represents  the  initial  indication 
that  nickel  in  this  type  of  experiment  would  cause  a non-zero  signal. 

1_ 


TABLE  D.1 .--Sunmary  of  Data  on  Shock  Demagnetization 


<^(mv)  , 4>{my) 


FT ^ 
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However,  it  was  thought  that  the  compression  of  the  teflon  insulator  might  ^ 

contribute  to  the  non-zero  signal  through  a dielectric  polarization  mechanism. 

76-037.  In  order  to  eliminate  the  ambiguity  with  respect  to  the 
I insulator  in  shot  76-032,  the  geometry  of  this  experiment  was  altered  as 

I ■ shown  in  Fig.  D.2.  The  return  path  for  the  signal  is  completely  independent 

t 

i of  the  shock  guard  ring  system,  including  the  insulator.  Direct  measurement 

' of  any  shock  polarization  is  thus  eliminated.  Again  a negative  non-steady 

emf  was  recorded  with  an  initial  value  of  zero.  (See  Fig.  D.l.)  The  time 
variation  of  the  signal  is  entirely  different  from  that  of  shot  76-032  for 
reasons  unknown. 

76-039.  This  experiment  was  designed  to  duplicate  76-037  except  for 
the  projectile  design.  The  projectile  used  was  machined  from  solid  nylon  in 
order  to  modify  the  pre-impact  noise  observed  in  this  type  of  experiment  in 
which  the  circuit  is  completed  at  impact.  The  resulting  voltage-time  profile, 
shown  in  Fig.  D.l,  is  again  initially  zero  but  positive  thereafter.  This 
change  in  sign  may  be  a result  of  the  lack  of  characterization  of  the  material 
with  respect  to  magnetic  properties. 

76-040.  The  guard  ring  was  completely  eliminated  in  this  experiment 
to  clearly  indicate  an  effect  involving  only  the  nickel  active  element  (as 
' seen  in  Fig.  D.3).  The  resulting  voltage-time  profile  is  again  positive 

(see  Fig.  D.l ). 

The  next  series  of  experiments  was  designed  to  show  specifically 
whether  or  not  the  observed  non-zero  signals  in  this  type  experiment  were 
due  to  demagnetization. 

76-051 . The  geometry  used  for  this  experiment,  where  the  nickel 


element  is  impacted  directly  by  an  aluminum  projectile,  is  shown  in  Fig.  D.4. 
This  shot  was  done  as  a control  experiment  for  the  shot  to  follow,  76-054, 
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Fig.  D. 2. --Modified  Concentric  Ring  Configuration  Used 
in  Shots  76-037  and  76-039. 
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Fig.  D. 3. —Experimental  Configuration  Used  in  Shot  76-040. 


and  the  modification  in  design  from  previous  shots  was  for  simplicity  of 
construction  and  analysis.  The  magnetization  state  was  not  altered  from  its 
supplied  configuration,  and  was  characterized  by  measuring  the  magnetic  field 
just  outside  the  body.  The  magnetic  induction  B was  approximately  0.8 
Gauss  with  components  in  directions  parallel  and  perpendicular  to  the  shock 
propagation  direction.  The  resulting  data  trace  has  much  the  same  character 
as  those  of  the  previous  series  with  an  amplitude  of  -40  mv.  (See  Fig.  0.3.) 

76-054.  In  this  experiment,  the  magnetization  of  the  nickel  was 
decreased  by  careful  demagnetization,  using  an  alternating  field,  to  a value 
less  than  0.05  Gauss  in  any  direction.  In  all  respects,  the  experiment  was 
conducted  in  the  same  manner  as  shot  76-051.  The  resulting  amplitude  varied 
between  3 mv  and  7 mv  with  an  initial  value  of  zero  (see  Fig.  0.5).  In  these 
two  experiments,  a reduction  in  initial  magnetization  resulted  in  a correspon- 
ding reduction  in  the  absolute  value  of  the  emf  recorded. 

76-062.  This  experiment  involved  the  use  of  a concentric  guard  ring 
system  as  used  in  shot  76-032  and  shown  in  Fig.  C.l.  The  material  used  for 
the  active  center  element  was  Nichrome  V,*  an  alloy  composed  of  80%  nickel 
and  20%  chromium.  The  material  was  chosen  because  of  the  mechanical  simi- 
larity to  nickel,  and  because  it  is  paramagnetic  at  room  temperature  and 
pressure.  The  voltage-time  profile,  shown  in  Fig.  0.5,  indicates  no  emf  pro- 
duced during  the  time  of  shock  transit  through  the  target  to  within  the 
recording  capability  of  the  instrumentation  system  of  0.01  mv. 

The  final  experiment  to  be  considered,  76-052,  was  designed  to  dupli- 

19 

cate  the  conditions  in  the  experiments  reported  by  Buzhinskii  and  Samylov. 

The  geometry  involved  is  the  same  as  that  used  to  investigate  the  effects  of 
surface  preparation,  and  is  discussed  in  Appendix  E (see  Fig.  E.l).  The  emf 

‘Trademark  of  Driver  Harris  Co.,  Harrison,  New  Jersey. 
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recorded  is  shown  in  Fig.  D.5.  It  does  not  appear  to  be  possible  at  this 
point  to  determine  how  the  signals  recorded  in  this  way  are  affected  by  the 
demagnetization.  This  determination  would  involve  separating  the  demagneti- 
zation effect  from  the  normal  response  of  a bimetallic  junction  under  these 
conditions. 
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APPENDIX  £ 

CALCULATIONS  AND  EXPERIMENTS  INVOLVING  THE 
SURFACE  LAYER  AT  THE  JUNCTION 


Consider  the  temperature  distribution  of  Fig.  3.4.  The  time  depen- 
dence of  the  initial  distribution  can  be  studied  using  the  one-dimensional 
heat  flow  equation.  The  appropriate  solutions  are  most  easily  found  using 
Green's  functions.  In  one  dimension,  the  diffusion  equation  for  temperature 
is 


9^T  1 aT 

3x2  - D3t 

where 


(E.l) 


< = thermal  conductivity 
p = mass  density 

C^  = constant  volume  specific  heat. 

The  unit  source  function  or  Green's  function  for  a unit  source  at  position 
x'  is  given  by  Carslaw  and  Jaeger^ ^ as 

G(D,  X.  X'.  t)  = — L_exp(-  . (E.2) 

In  a situation  where  two  different  metals  with  associated  thermal 
properties  are  involved,  as  in  Fig.  3.4,  a more  complex  source  function  must 
be  used.  Carslaw  and  Jaeger^ ^ give  the  Green's  function  for  a source  in  a 


region  2 located  a distance  H from  the  interface  at  the  origin  as 
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where 

and 

3 = 1/a 

The  effect  of  the  interface  on  a source  in  region  1 is  that  in  region 
2 the  source  appears  displaced  to  a position  3Jl  and  reduced  in  amplitude  by 
a factor  23/ (1+6),  while  in  region  1 an  additional  source  appears  in  a 
reflected  location  and  of  amplitude  (1-3/1+8).  A similar  situation  exists 
for  a source  in  region  2.  Now  to  find  the  temperature  history  due  to  a square 
initial  temperature  profile  we  have: 


T(x.t) 


0 

Tg  ’^G{D2,  X,  3)1.  t)d)t 

-a 
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a 

+ jlgCGCDg,  X,  t)  + G(D2,  X,  -I,  t)]d)l  for  'x  > 0 . (E.5) 

0 


For  simplicity  is  taken  to  be  the  zero  of  temperature.  A similar  expres- 
sion holds  for  X < 0;  however,  to  find  T(x  = 0,  t)  one  need  only  solve  for 
one  half  space.  The  solution  is  given  in  terms  of  error  functions  as 


T(x,t)  = 


T f 
'2  ^ 1+3 


erf  ( 


x+Sa 


;)  - y erf  ( 


x-a  y 


+ J > for  X > 0 . (E.6) 

^ “ AD^ 

Of  primary  interest  is  the  solution  at  x = 0 which  is 

T(O.t)  = erf  (-^)  + ^ erf  (-^) } (E.7) 

If  region  1 is  taken  to  be  constantan  and  region  2 is  copper,  T(0,t)  becomes 
T(0,t)  = T2{^erf  (2/F7t)  + | erf  /?7t)}  (E.8) 

2 

where  x = a /D2. 

Two  experiments  were  conducted  to  determine  the  effects  of  surface 
preparation  in  a geometry  in  which  the  two  surfaces  making  up  the  junction 
are  held  in  contact  before  shock  compression.  The  configuration  for  each  of 
these  experiments  is  given  in  Fig.  E.l.  In  these  two  experiments  the  only 
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controlled  difference  was  the  surface  on  both  the  copper  buffer  and  the  con- 
stantan  element.  The  materials  used  for  the  junction  in  both  shots  were  B- 
constantan  and  OFHC  copper  as  described  in  section  4.2.  In  shot  77-038,  a 
highly  polished  surface  was  prepared  using  a series  of  diamond  polishing 
steps  concluded  with  a final  polish  using  0.05  y alumina.  Attention  was  given 
throughout  to  maintain  flatness,  and  this  was  achieved  to  within  ten  wave- 
lengths of  light  over  90%  of  the  area.  In  shot  77-037,  lapped  surfaces  were 
prepared  using  a Lapmaster*  lap  with  #1700  lapping  compound.  This  surface  is 
believed  to  have  damage  on  the  order  of  10  y and  perhaps  as  much  as  30  y in 

depth.  The  results  of  these  two  experiments,  as  well  as  those  reported  by 
28 

Nesterenko,  are  included  in  Table  E.l.  The  "measured  temperature"  listed 
in  the  table  is  to  be  interpreted  as  the  temperature  that  would  be  needed  for 
a standard  pressure  thermocouple  of  the  respective  materials  to  produce  the 
measured  emf.  The  voltage-time  profiles  for  shots  77-038  and  77-039  are 
given  in  Fig.  E.2. 
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TABLE  E.l .--Summary  of  Data  in  Surface  Preparation  Studies 


^(mv) 


26 


t(/isec) 


0.0  0.5  1.0 


Fig.  E.2.— Voltage-Time  Profil 


t(^sec) 

0.0  0.5  1.0 


for  Shots  77-038  and  77-039. 
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APPENDIX  F 

THE  EFFECTS  OF  A CHANGING  INDUCTANCE  AND  SPECIAL  RELATIVITY 
WHEN  RADIAL  PRESSURE  RELIEF  IS  ALLOWED 

Consider  the  effect  of  the  changing  inductance  in  a thermoelectric 
circuit  as  a possible  source  of  the  ramping  discussed  in  section  5.1.  The 
pressure  relief  waves  behind  the  shock  will  accelerate  particles  radially 
outward.  The  radius  of  the  constantan  will  thus  be  increasing  with  time.  As 
a model  to  investigate  this  problem  consider  a circuit  with  a constant  emf 
source,  ((>,  representing  the  thermal  emf,  a resistive  element,  R,  representing 
the  input  impedance  of  the  recording  circuit,  and  an  inductance,  L.  The 
inductance  will  have  a component,  L.| , which  changes  with  time  representing 
the  changing  inductance  of  the  expanding  conductor.  The  differential  equation 
for  the  circuit  is 

HT  ‘IL, 

I-  f ’ -ar* 

At  t ==  0 we  will  assume  the  current  has  a value  (p/R.  This  equation  has 
the  solution 


where 


Ae 


(-t/x) 


R + 


dt 


dL, 


■ar 


R^  + 


dt 


(F.2) 


1 


A * 


and  T 


dL, 
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To  estimate  the  value  of  we  will  model  the  expanding  constantan 

cylinder  by  the  solid  center  conductor  of  a coaxial  transmission  line.  The 
inductance  of  a coaxial  conductor  of  length  in  MKS  units  is 


2, 


(F.3) 


where  is  the  permeability  in  a vacuum,  y is  the  permeability  of  the 
conductor  material,  a is  the  radius  of  the  center  conductor  and  b is  the 
radius  of  the  outer  conductor.  Now  if  a increases  with  time,  but  b 
remains  constant,  the  rate  of  change  of  the  inductance  is 


dLi 

dt 


V 1 ^ 

2Tr  a dt 


(F.4) 


If  a and  I are  the  radius  and  length  (10  mm,  and  5 mm)  of  the  constantan 

and  ^ has  a value  of  1 mm/ysec  then  has  the  value  1 x 10'^  henry/sec. 

Since  Eq.  (F.3)  is  derived  assuming  an  infinitely  long  coaxial  transmission 

line  the  value  of  L for  the  dimensions  here  will  be  inaccurate.  However, 

dL, 

this  calculation  should  give  a reasonably  good  approximation  for  . 

The  time  constant  x is  0.5  ns  where  L = 0.05  yhenry,  R = 100  n, 

^^1  -4 

and  = 1 X 10  henry/sec.  The  value  for  R is  the  actual  input  impe- 
dance for  our  recording  system  and  the  value  for  L is  an  estimate  of  the 
total  inductance  of  the  experimental  configuration.  The  effect  on  t of  the 
changing  inductance  is  negligible. 

The  measured  steady  state  emf  will  be 


m 


1 

1 4.  J L 

' R dt 


Using  a binomial  expansion  this  becomes 


129 


T 


! I 


dL 


(,  .1^+  (ldk)2  * 

R dt 


dt^ 


. . ) 


(F.5) 


1 


The  first  term  in  has  a value  of  1 x 10“®  for  the  values  given  above 

which  is  obviously  negligible  compared  to  1. 

Although  the  assumptions  made  here  to  do  this  calculation  are  not 

dL, 

valid,  the  estimate  of  the  value  of  should  be  within  1 or  2 orders  of 

magnitude.  The  effect  of  this  changing  inductance  would  still  be  negligible 
and  could  not  account  for  the  observed  ramping. 

To  examine  the  relativistic  effect  on  the  electric  field  due  to 
radial  pressure  relief  in  the  experiments  discussed  in  section  5.1  consider 
a conducting  cylinder  with  a uniform  current  density  I,  and  electric  field 
E in  the  direction  of  the  cylinder  axis.  The  pressure  relief  waves  from 
the  edges  will  accelerate  the  particles  radially  so  that  the  velocity  v will 
oe  perpendicular  to  the  electric  field  everywhere.  If  the  electric  field  is 
locally  unchanged  in  the  moving  material,  the  electric  field  measured  in  the 
stationary  laboratory  frame  is 


E = y(E'  - B') 


where  E'  and  B'  are  measured  in  the  moving  frame,  and 

1 


(F.6) 


(F.7) 


For  a uniform  radial  particle  velocity  of  1 mm/ysec,  which  is  larger  than 
would  be  possible  in  any  experiments  done  here,  the  value  of  y is 
1.00000167.  The  term  in  B'  is  negligible  compared  to  E'  so  that  E and 
E'  are  indistinguishable  in  this  type  of  experiment. 


